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ABSTRACT 

An  experimental  investigation  was  conducted  of  crater  growth 
and  ejecta  cloud  formation  from  the  impact  of  3.  Z  mm  aluminum 
spheres  on  thick  aluminum  targets  at  7  km/sec.  Crater  growth  and 
transient  shape  were  determined  through  sequential  flash  x-rays. 
Growth  followed  a  decaying  exponential  pattern,  and  the.  tangent  angle 
to  the  crater  wall  at  the  target  surface  remained  virtually  constant. 
Relationships  between  ejecta-cloud  parameters  and  crater  diameters 
were  investigated.  Cloud-edge  motion  was  determined  and  an  effort 
made  to  determine  particle  origin.  Velocities  of  discrete  particles 
in  the  cloud  were  determined.  No  direct  relationship  between  cloud 
parameters  and  crater  dimensions  could  be  established. 
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INVESTIGATION  OF  CRATER.  GROWTH  AND 
EJECTA  CLOUD  RESULTING  FROM 
HYPERVELOCITY  IMPACT  OF  ALUMINUM 
SPHE.UES  ON  THICK  ALUMINUM  TARGETS 


I.  Introduction 


Background 

With  the  advent  of  the  space  age,  the  problem  of  material 
damage  caused  by  projectiles  impacting  targets  at  hypervelocities 
became  one  of  paramount  concern.  Earlier  calculations,  based  on 
Bernoulli's  equation  for  incompressible  fluid  flow,  were  carried  out 
by  Pugh  in  the  United  States  and  by  Hill  et  M.  in  Great  Britain  in 
connection  with  the  "shaped-charge"  developed  in  World  War  II  for 
penetrating  armor.  These  calculations  were  found  to  be  inapplicable 
to  the  problem  of  randomly  oriented  bodies  impacting  solid  targets. 

In  1968,  employing  classical  hydrodynamics  of  compressible  media 
as  a  starting  point,  Bjork  developed  what  is  generally  considered  the 
first  comprehensive  theoretical  treatment  of  the  problem  of  hyper- 
velocity  impact  and  resulting  crater  formation.  Bjork  used  a  digital 
computer  to  obtain  numerical  solutions  to  the  two-dimtnsional,  time- 
dependent  hydrodynamic  equations  for  a  compressible  medium 
(Ref  1;220).  Subsequently,  the  hydrodynamic  approach  was  employed 
by  Walsh  and  Tillotson  (Ref  30)  and  by  Riney  (Ref  24). 
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Concurrently  with  the  development  of  these  theories,  experi¬ 
menters  in  laboratories  both  here  and  abroad  started  amassing  data 
which  might  presumably  serve  to  support  or  refute  these  or  other 
theories.  An  idea  of  the  acceleration  of  the  pace  of  theory  develop¬ 
ment  and  experimentation  in  hypervelocity  impact  can  be  inferred  from 
the  volumes  reporting  the  proceedings  of  the  biennial  Hypervelocity 
Impact  Symposia.  The  unclassified  results  of  the  Second  Symposium 
occupy  a  single  volume,  while  the  unclassified  proceedings  of  the 
Seventh  Symposium  fill  six  separate  volumes,  each  of  size  comparable 
to  that  of  the  single  volume  resulting  from  the  Second  Symposium. 

Unfortunately,  much  of  the  data  thus  acquired  and  presented  is 

not  suitable  for  completely  testing  the  cited  theories.  The  practical 

limit  on  velocities  attainable  in  light-gas  guns  launching  projectiles  of 

known  dimensions  and  mass  is  currently  on  the  order  of  10  km/sec. 

Such  impacts  will  generate  and  sustain  pressures  sufficient  to  ensure 

true  hydrodynamic  conditions  for  only  a  few  microseconds.  Thus  a 

large  part  of  the  crater  formation  process  in  laboratory  experiments 

occurs  at  lower  pressures  where  material  properties  cannot  be 

ignored.  While  conceding  that  in  the  final  stages  of  crater  formation 

material  characteristics  may  become  significant,  Olshaker  and  Bjork 

dismiss  the  problem  of  determining  precisely  at  what  point  in  the 

process  this  transition  occurs  in  the  belief  that  at  the  higher  meteoroid 

velocities  of  interest  (11  to  72  km/sec),  hydrodynamic  effects  will  far 

outweigh  any  late  stage  strength  effects  (Ref  19:225).  The  treatment 
l 


2 


of  the  problem  by  Walsh  and  by  Riney  is  different  in  that  both  ter¬ 
minate  the  hydrodynamic  computations  well  before  the  crater  has 
stopped  growing.  Implicit  in  both  treatments  is  the  assumption  that 
the  strength  of  the  impacted  material  cannot  be  ignored,  even  at 
arbitrarily  high  impacting  velocities  (Ref  11:239). 

From  their  observations,  most  experimenters  echo  Riney' s 
and  Walsh's  insistence  on  including  consideration  of  target  strength 
in  any  tenable  theory  of  crater  formation.  Rae  and  Kirchner  state 
the  majority  view  succinctly  when  they  write: 

.  .  .  The  establishment  of  a  crater  of  fixed  size  im¬ 
plies  that  material  has  been  brought  to  rest,  and.  .  . 
there  is  no  mechanism  for  accomplishing  this  feat 
within  the  framework  of  an  inviscid  theory.  Thus  it 
appears  that  at  large  time  a  transition  must  be  made  to  a 
theory  which  accounts  for  the  strength  of  the  target.  .  .  . 

Thus  we  ought  to  assign  as  a  boundary  of  the  hydro- 
dynamic  theory  some  level  of  pressure  comparable 
with  target  strength  (Ref  23:210). 

Precisely  what  this  boundary  represents  is,  unfortunately,  not 
at  all  clear.  That  it  is  in  some  way  related  to  target  strength  is  quite 
well  established.  Several  investigators  have  found  that  cratering 
efficiency  in  terms  of  projectile  energy  per  unit  crater  volume  is 
directly  proportional  to  the  Brinell  hardness  number  of  the  target 
material  over  a  wide  range  of  materials  and  values  of  hardness 
(Ref  10:165).  The  influence  of  the  density  of  either  the  target  or  the 
projectile  remains  a  subject  of  considerable  controversy.  Expressed 
opinions  range  from  quadratic  relations  (between  crater  volume  and 
projectile  density)  to  no  effect  (Ref  4:334).  Frasier  et  find  good 
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agreement  between  Walsh's  hydrodynamic  treatment  and  their  experi¬ 
mental  data  down  to  pressures  about  an  order  of  magnitude  above  the 
Hugoniot  elastic  limit  of  the  target  (Ref  8).  Sorensen  concludes  that 
all  metal-on-metal  impacts  up  to  7.  5  km/sec  produce  crater  volumes 
which  are  inversely  proportional  to  target  shear  strength  to  the  0.  845 
power  (Ref  26).  Several  investigators  employ  the  ratio  of  projectile 
velocity  to  target  dilational  wave  velocity  (v/c)  in  the  empirical 
equations  which  they  derive,  but  there  is  little  agreement  as  to  how 

this  ratio  enters  the  equation.  Maiden  find  that  penetration 

2 

varies  as  (v/c)  (Ref  18);  on  the  other  hand.  Summers  £t  al.  conclude 
that  it  varies  as  (  v/c)^^^  (Ref  27). 

On  the  fringes  of  the  empiricist  group  stand  a  few  who  deny  the 
validity  of  the  hydrodynamic  approach  in  toto.  For  example,  Feldman 
finds  that  the  ratio  of  kinetic  energy  to  crater  volume  is  the  same  for 
jets  and  pellets  and  therefore  concludes  that  hypervelocity  impact  is 
basically  a  problem  In  plastic  flow,  not  shock  wave  phenomena  (Ref  7). 
Engle  proposes  that  hypervelocity  impact  is  best  treated  by  the  analogy 
of  liquid  droplets  impacting  a  fluid  surface  (Ref  6).  And  Pond  et  al. 
find  that  major  proportions  of  kinetic  energy  are  expended  in  strain- 
hardening  various  metals,  a  viewpoint  which  is  in  marked  contrast  to 
the  hydrodynamic  approach  (Ref  22). 

But  by  far  the  majority  of  experimenters  share  the  more 
moderate  view  that  although  hypervelocity  impact  may  be  amenable  to 
hydrodynamic  treatment,  still  for  impact  velocities  at  least  up  to  the 
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limits  currently  attainable,  target  strength  is  significant  in  determin¬ 
ing  final  crater  dimensions.  In  fact,  Kineke  and  Richards  find  no 
indication  that  the  influence  of  mechanical  strength  is  decreasing  for 
impacts  up  to  15.  5  km/sec  (Ref  15). 

A  report  published  by  Hermann  andJones  in  an  effort  to  analyze 
and  correlate  the  various  experimental  data  aptly  summarizes  the 
state  of  confusion  currently  existing: 

A  great  deal  of  experimental  data  has  been  published 
by  a  number  of  laboratories,  and  considerable  theoret¬ 
ical  work  has  appeared  during  the  last  few  years. 

Limited  comparisons  of  data  from  one  laboratory  with 
those  from  another  and  of  experimental  data  with 
theoretical  predictions  have  shown  some  large  dis¬ 
crepancies.  Considerable  confusion  exists  because, 
except  in  a  few  instances,  the  data  from  one  laboratory 
are  not  directly  comparable  with  those  from  another 
since  identical  materials,  projectile  shapes,  and 
velocity  ranges  were  not  used.  Each  laboratory  has 
produced  a  different  empirical  expression  which  was 
found  to  fit  the  limited  range  in  experimental  parameters 
explored  at  the  laboratory.  The  empirical  expressions 
are  more  or  less  contradictory,  and  when  extrapolated 
to  velocities  of  interest  in  space  applications,  lead  to 
large  disagreements  in  predicted  penetrations 
(Ref  12:390), 

Although  this  particular  report  was  published  some  five  years  ago,  the 
situation  does  not  appear  to  have  improved  significantly  since  then. 

In  his  Summary  of  Theoretical  and  Experimental  Studies  of 
Crater  Formation  presented  at  the  Sixth  Hypervelocity  Impact  Sym¬ 
posium,  Eichelberger  makes  some  penetrating  (and  occasionally 
caustic)  comments  on  this  apparent  inability  of  theoreticians  and 
experimentalists  to  find  a  com.mon  meeting  ground  to  permit  the 
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testing  of  theory  with  meaningful  experimental  data.  He  concedes  that 
because  of  velocity  limitations,  the  early  or  purely  hydrodynamic 
stages  will  probably  not  vary  sufficiently  to  permit  a  clearcut  choice 
between  theories.  Rather,  he  says,  "It  is  more  likely  that  experi¬ 
ments  involving  detailed  observations  of  transient  conditions  during 
the  later  stages  of  the  crater  formation  process  will  provide  decisive 
comparisons"  (Ref  5:688-689). 

Purpose 

In  all  the  welter  of  reported  observational  data,  surprisingly 
little  attention  has  been  paid  to  the  cloud  of  ejecta  which  inevitably 
accompanies  and  results  from  the  process  of  crater  formation. 
Apparently  most  experimenters  have  considered  the  cloud  only  an 
annoying  obstacle  preventing  their  direct  viewing  of  the  crater  forma¬ 
tion  process  in  which  they  were  primarily  interested.  Yet  it  seems 
plausible  that  since  the  ejecta  cloud  is  intimately  associated  with 
whatever  process  or  processes  are  operative  in  forming  the  crater, 
valuable  information  about  these  processes  might  be  obtained  by 
careful  examination  and  analysis  of  pictures  taken  of  the  ejecta  cloud 
during  the  cratering  process. 

If,  as  Riney  suggests,  the  mechanism  of  crater  formation  is 
essentially  one  of  cavitation  (Ref  24;lbl),  then  determination  of  the 
trajectories  of  the  particles  thus  ejected  should  enable  one  to  acquire 
some  information  on  the  forces  and  conditions  responsible  for  their 
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ejection.  Kinalow's  experimental  observations  support  Riney's  views. 
He  infers  that  as  the  expanding  shock  wave  becomes  detached  from  the 
crater,  material  flows  along  the  walls  of  the  crater  and  is  ejected  at 
velocities  up  to  eight  or  ten  times  the  impacting  velocity  of  the  pro¬ 
jectile  (Ref  16:279).  While  agreeing  in  principle  with  this  view,  most 
other  experimenters  place  a  lower  value  on  the  maximum  eiecta 
velocity  of  perhaps  triple  the  impact  velocity  (Refs  14;  17).  Klneke 
also  notes  that  as  ejection  continues,  the  '’ected  pieces  increase  in 
size  and  decrease  in  velocity  (Ref  14:351-  352).  Detailed  studies  of  the 
ejecta  produced  by  hypervelocity  impact  of  metallic  spheres  into  rock 
targets  conducted  by  Gault  et  al.  confirm  this  finding  of  generally 
monotoni^ally  increasing  ejecta  size  combined  with  decreasing  ejecta 
velocities.  Gault  also  finds  that  for  basalt  ejecta,  following  the  initial 
jetting  phase,  the  angle  of  ejection  (measxired  from  the  target  face) 
increases  to  60  degrees,  decreases  to  about  50  degrees,  and  finally 
tends  toward  90  degrees  (Ref  9:446-449).  This  departure  from  a 
rnonotonicaily  increasing  ejection  angle  has  apparently  not  been 
observed  or  reported  in  the  ejecta  from  metallic  targets.  But  as 
noted  previousiy,  the  ejecta  patterns  from  metallic  targets  have  not 
been  studied  as  extensively  as  has  basalt  ejecta. 

These  observed  variations  in  ejecta  parameters  with  time  and 
their  inferred  dependence  on  conditions  existing  within  the  crater  at 
the  time  of  ejection  reinforce  the  view  that  a  detailed  examination  and 
analysis  of  sequential  pictures  of  the  ejecta  cloud  might  provide 
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Information  on  crater  parameters  as  a  function  of  time  during  the 
crater  formation  process.  The  purpose  of  this  study,  then,  was  to 
examine  possible  relationships  between  ejecta  cloud  parameters  and 
crater  growth  and  shape. 

The  study  was  designed  in  two  independent  but  closely  related 
phases.  In  a  series  of  hypervelocity  impact  experiments  crater  growth 
data  would  be  obtained  through  the  use  of  a  battery  of  sequential  flash 
x-rays.  At  the  same  time,  optical  photographs  would  be  made  of  the 
ejecta  cloud  resulting  from  the  impact.  It  was  felt  that  accurate 
measurement  of  changing  cloud-edge  and  ejecta  particle  positions  in 
sequential  optical  photographs  would  yield  information  on  ejecta 
particle  velocities  and  times  and  points  of  ejection.  Concurrently, 
sequential  x-ray  photographs  (radiographs)  of  crater  profile  and  plan 
would  provide  information  on  crater  parameters  and  shapes  as  a 
function  of  time.  Comparison  of  ejecta  parameters  determined  from 
cloud  photographs  with  radiographic  data  on  crater  growth  would 
permit  correlation  of  the  two  phenomena.  The  establishment  of  a 
significant  correlation  would  provide  an  extremely  useful  technique  to 
experimentalists  concerned  with  the  problem  of  testing  theoretical 
predictions  of  crater  growth  in  thick  metallic  targets. 
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II.  Postulated  Model  and  Experimental  Approach 

Postulated  Model  for  CratertnR  in  Thick  Targets 

In  order  to  establish  logical  relationships  between  crater 
formation  and  the  resulting  ejecta  cloud,  it  was  necessary  to  postulate 
a  model  of  crater  formation.  Despite  the  cited  disagreements  on  the 
relative  importance  and  effects  of  material  properties,  most  investi¬ 
gators  agree  that  crater  formation  in  thick  targets  occurs  in  the 
following  stages.  The  projectile  penetrates  the  target  surface  crea¬ 
ting  a  very  intense  shock  wave.  Cavitation  is  initiated  behind  the 
shock  wave.  The  velocity  of  the  expanding  crater  surface  decreases, 
and  the  shock  wave  detaches  from  the  surface.  Material  flows  along 
the  crater  walla  and  is  ejected.  The  crater  continues  to  expand  av  a 
decreasing  rate  until  crater  growth  is  arrested  by  the  dynamic 
strength  of  the  material.  The  amplitude  of  the  expanding  shock  wave 
decreases  as  energy  is  dissipated  throughout  the  material  until  the 
wave  decays  into  an  elastic  wave  (Ref  16;Z79). 

A  simplified  schematic  of  the  outlined  cratering  process  is 
presented  in  Fig.  1.  Consideration  of  the  previous  discussion  and 
crater  geometry  leads  to  the  following  set  of  assumptions  constituting 
a  proposed  model  for  crater  and  ejecta  cloud  formation; 

a.  Material  departs  the  crater  as  ejecta  along  vectors  tangent 
to  the  crater  wall  at  or  very  close  to  the  original  target 
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Fig.  1.  Schematic  of  Hyperv,;locity  Impact  Cratering  Procesa 

surface.  Therefore,  the  points  of  particle  ejection  are 
Virtually  identical  to  the  crater  radius  at  all  tirues, 

b.  At  any  instant,  particles  are  being  ejected  at  nearly 
identical  velocities. 

c.  As  the  crater  enlarges,  tlie  angle  of  ejection  increases 
toward  the  vertical  with  respect  to  the  target  face. 

d.  As  the  pressure  dissipates,  the  velocity  of  ejection 
decreases. 

e.  No  significant  pressure  field  exists  outside  the  crater  to 
alter  particle  velocity.  Therefore,  after  being  ejected, 
particles  travel  in  straight  lines  at  constant  velocities 
except  for  deceleration  due  to  atmospheric  crag  effects. 


10 


Approach  to  Testing  the  Model 


As  the  ejection  process  occurs,  the  characteristic  shape  of 
the  ejecta  cloud  varies  with  time.  For  purposes  of  the  following 

discussion,  "early  time"  refers  to  times  less  than  approximately 

/ 

10-15  microseconds  after  impact.  The  remainder  of  the  ejection 
process  will  be  referred  to  as  "late  time".  In  early  time,  the  cloud 
edges  form  smooth  curves  which  slope  inward  toward  the  cloud 
centerline  from  their  intersection  with  the  target  face.  After  reaching 
a  limiting  point  of  minimum  width  of  the  "throat"  of  the  cloud,  the 
cloud  edges  slope  back  away  from  the  cloud  centerline  and  become 
less  distinct  with  Increasing  distance  from  the  target  face  (Fig.  2). 

To  a  first  approximation,  the  cloud  appears  symmetrical  about  a 
centerline  which  presumably  coincides  with  the  trajectory  of  the 
impacting  projectile.  At  this  stage  in  the  process,  discrete  particles 
forming  the  ejecta  cloud  are  not  discernible.  However,  it  would  seem 
logical  that  the  cloud  edges  constitute  lines  of  discrete  although 
indistinguishable  particles.  If  the  velocities  of  ejection  are  monoton- 
ically  decreasing  at  this  stage,  and  if  the  particles  can  be  assumed  to 
travel  at  relatively  constant  velocity  following  their  ejection,  then  it 
follows  that  the  same  particles  which  form  the  cloud  edge  at  one  point 
in  time  must  necessarily  form  the  cloud  edge  in  subsequent  times. 
Hence  it  should  be  possible  to  use  a  set  of  sequenced  photographs  of  a 
developing  ejecta  plume  to  reconstruct  the  straight-line  trajectories 
of  individual  particles.  The  distance  traveled  by  a  particle  between 


11 


Fig,  2.  "^pical  Early-Time  Photograph  of  Ejecta  Cloud 
^/Magnification  =5.3  (nom);  Time  =  1.03^  sec; 

Velocity  of  Impact  =  7  'Kni/se£J 

cloud  photographs  is  in  the  same  ratio  as  the  time  between  frames. 
Once  the  correct  geometrical  relationships  are  established,  particle 
velocities  can  be  uniquely  deterniined.  By  extending  the  straight  line 
of  each  particle  trajectory  back  to  tne  target  face,  its  origin,  ejection 
time,  and  angle  of  ejection  can  all  be  determined. 

In  order  to  test  the  feasibility  of  establishing  the  required 
geometrical  relationships,  a  simpjc  computer  program  was  written  to 
eject  hypothetical  particles  which  travel  at  constant  velocity.  Particle 
origin  was  moved  outward  with  time  (siinulating  crater  growth). 
Velocities  of  ejection  were  monotonicalJ y  decreased,  and  angles  of 
ejection  from  the  target  face  were  monotonically  increased.  A  plot  of 
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particle  positions  thus  determined  as  a  function  of  time  for  three  times 
appears  as  Fig.  3,  with  straight  lines  indicating  particle  trajectories. 

The  problem  of  constructing  the  particle-trajectory  straight 
lines  by  working  backward  from  cloud  photographs  was  Initially 
attempted  graphically.  However,  as  can  be  seen  from  Fig.  3,  the 
trajectory  lines  intersect  the  cloud  edge  curves  at  such  low  angles  that 
a  graphical  solution  is  subject  to  considerable  error.  Consequently,  a 
computer  program,  was  developed  to  solve  the  problem  by  repeated 
iteration. 

The  cloud  edges  were  fit  by  a  standard  least  squares  program 
to  polynomials  of  order  2  through  10,  and  the  best  fit  selected  to 
represent  the  cloud  edge  in  the  program.  The  x-axis  was  placed 


Fig.  3.  Simulated  Ejecta  Curves  Produced  by  Constant- 
Velocity  Particle  Ejection  Program 


13 


along  the  symmetry  axis  of  the  cloud,  with  x-distance  increasing  away 
from  the  target  face  (uprange).  The  y-axis  was  placed  on  the  target 
face,  with  origin  at  the  intersection  of  the  x-axis  with  the  target  face. 
With  the  three  cloud  curves  designated  A,  B,  and  C  in  order  of 
increasing  time,  a  point  was  arbitrarily  selected  on  the  A  curve.  A 
logical  estimate  was  made  of  the  position  of  origin  on  the  target  face. 
The  point  of  intersection  with  the  B  curve  of  a  straight  line  from  the 
assumed  origin  through  the  point  on  the  A  curve  was  then  determined 
by  standard  Newton-Raphson  iteration  techniques.  The  Intersection 
of  this  line  with  the  C  curve  was  then  similarly  determined.  The 
distances  between  curve  intersections  were  compared  with  the  time 
differences  between  the  frames  represented  by  the  three  curves.  On 
the  basis  of  the  resultant  error,  the  assumed  origin  was  corrected  to 
bring  the  distance  ratios  into  agreement  witli  the  interframe  times, 
and  the  iteration  process  continued  until  acceptable  agreement  was 
achieved.  A  graphical  portrayal  of  the  scheme  is  presented  in  Fig.  4. 

The  accuracy  of  the  completed  program  (Appendix  A)  was 
checked  against  the  artificial  cloud  generation  program  with  the 
results  shown  in  Table  I.  Since  the  resulting  errors  appeared  to  be 
within  reasonable  limits  and  considerably  less  than  might  be  expected 
from  a  purely  graphical  solution,  use  of  the  program  was  considered 
warranted. 

In  later  times  during  the  ejection  process,  as  increasingly 
larger  particles  are  ejected,  the  characteristic  ejecta  cloud  shape 
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Fig.  4.  Graphical  Presentation  of  Cloud-Edge 
Computer  Program  Routine 


undergoes  marked  alterations  (Fig.  5),  The  bottom,  inward- sloping 
segment  of  the  cloud  edge  gradually  disappears  until  the  narrowest 
portion  or  throat  of  the  cloud  coincides  with  its  intersection  with  the 
target  face.  The  cloud  now  consists  largely  of  discrete  and  readily 
discernible  particles.  Inclusion  of  these  particles  in  the  formerly 
smooth  cloud  edge  causes  the  cloud  edge  to  appear  more  ragged. 

Rather  than  a  line  of  indistinguishable  particles  it  now  represents  a 
train  of  photographically  resolved  fragments.  Because  the  cloud  edge 
is  not  clearly  defined  and  because  its  shape  and  dimensions  vary  only 
a  slight  amount  in  the  later  stages  of  ejection,  it  was  not  felt  advisable 
to  attempt  to  determine  particle  trajectory  from  the  cloud  edges  in 
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TABLE  I 
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Fig.  5.  '^pical  Late-Time  Photograph  of  Ejecta  Cloud 

^Magnification  =  0.  34  (nom);  Time  =  37.  84  p.  secj 
Velocity  of  Impact  =  7  km/se£7 

this  portion  of  the  ejection  process.  Rather,  it  was  decided  to  attempt 
to  identify  and  follow  specific  particles  from  frame  to  frame. 

True  lateral  particle  motion  in  the  y-direction  is  recorded 
only  for  those  particles  on  the  cloud  edge.  (All  other  particles 
possess  a  velocity  vector  toward  or  away  from  the  camera.)  First 
consideration  would  then  indicate  that  only  particles  lying  on  the  cloud 
edge  could  be  used  for  velocity  determination.  However,  if  the 
particles  are  ejected  from  the  crater  wall  with  some  outward  velocity 
vector,  then  the  ejecta  cloud  must  at  all  times  be  hollow.  Experi¬ 
mental  evidence  supports  this  conclusion.  Ejecta  impacts  recorded 
on  the  uprange  wall  of  the  target  tank  invariably  reveal  a  roughly  j 
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circular  pattern,  with  virtually  no  evidence  of  significant  impacts 
within  this  circle.  This  pattern  of  ejecta  impact  was  recorded  on  a 
clean  piece  of  paper  affixed  to  the  downrange  face  of  the  vertical  x-ray 
cassette  (to  be  described  later)  (see  Fig.  6). 

If  the  cloud  is  indeed  hollow,  as  this  evidence  suggests,  then 
all  particles  visible  in  the  ejecta  cloud  must  in  fact  be  on  or  very  close 
to  the  edge  of  the  cloud.  And  if  symmetry  can  be  assumed,  the  radial 
motion  of  any  particle  between  frames  must  be  represented  by  the 
y-motion  of  the  cloud  edge  between  corresponding  x-posltions  in  the 
two  frames  (see  Fig.  7).  Under  these  assumptions,  the  true  velocity 
of  any  particle  can  be  determined,  and  this  approach  was  decided  upon 
for  late-time  frames. 
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Fig.  6,  Ejecta  Impact  Pattern  on  Paper  Placed  on  Face  of  Vertical 
X-ray  Cassette  Located  21.6  cm  Uprange  From  Target 


Fig.  7.  Graphical  Presentation  of  Scheme  for  Solution 
From  Discrete  Particle  Travel 
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Fig.  8.  Schematic  of  AFML  I.ight-Gas  Gun  Range 

The  flash  x-ray  systems  employed  in  the  range  facility  con¬ 
sisted  of  three  300  kv  and  two  105  kv  heads.  When  pulsed,  these 
heads  emit  a  very  intense  burst  of  x-ray  energy  of  approximately 
20-30  nanosecond  duration.  This  short  pulse  duration  produces 
radiographs  of  high-speed  events  with  negligible  motion  blur.  The 
output  of  x-ray  energy  from  each  head  was  sensed  by  a  switch 
physic^illy  mounted  on  each  head  (Ref  29).  The  electrical  pulse 
generated  by  each  switch  was  used  to  fire  a  xenon  flash  tube  (winker) 
which  was  photographed  on  the  range  time- recording  films  (to  be 
discussed  later).  This  pennitted  very  accurate  determination  of  the 
time  of  firing  of  each  x-ray  head.  Inasmuch  as  a  complete  description 
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of  the  facility  is  available  elsewhere  (Ref  28),  only  those  modifications 
of  and  additions  to  range  facilities  and  apparatus  necessary  for  this 
experiment  will  be  discussed  in  detail.  Radiography  and  optical 
photography  layouts  are  discussed  in  Chapters  III  and  IV. 

Time  Determination 

Elapsed  times  between  electronic  events  occurring  in  the 
course  of  hypervelocity  impact  experiments  conducted  on  the  AFML 
range  are  recorded  on  and  determined  from  reel-type  streak  camera 
(Fastax)*  filrris.  The  .  ^cording  of  the  time  of  non-electronic  events 
such  as  impact  is  more  difficult.  This  problem  was  solved  by  aiming 
a  photomultiplier  (PM)  tube  at  the  target  face.  The  impact  flash 
caused  the  PM  to  put  out  a  low-voltage  signal  which,  after  amplifica¬ 
tion,  was  similarly  used  to  fire  xenon  winkers. 

With  this  basic  system,  the  relative  time  between  any  two 
events  can  be  determined  to  an  accuracy  of  60  nanoseconds  (Ref  28:19). 
At  the  outset  of  the  experiment  this  basic  Fastax  winker  system  was 
used  to  determine  the  times  of  x-ray  firings,  the  output  of  each  x-ray 
head  switch  being  recorded  by  a  separate  winker.  Under  the  basic 
system,  these  times  could  be  related  to  the  time  of  impact  only  by 
employing  projectile  velocity.  In  an  effort  to  decrease  the  possible 
error  in  time  of  impact  resulting  from  either  errors  in  velocity  deter¬ 
mination  or  deceleration  of  the  projectile  after  such  determination,  an 

*  Wollensak  16  mm  "Fastax"  oscillographic  camera. 
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additional  rotating-drum  type  streak  g  camera  (Dynafax)*  was 
employed  to  record  the  actual  impact  flash  as  well  as  such  events  as 
x-ray  firings  and,  later,  PM  tube  output  and  camera  shutter  signals. 
The  Dynafax  is  capable  of  attaining  about  three  times  the  film  speed  of 
the  Fastax  camera.  By  essentially  "paralleling"  the  two  camera 
records,  it  was  possible  to  determine  relative  time  between  events 
with  a  maximum  error  of  about  40  nanoseconds. 

The  procedure  employed  to  determine  the  times  of  image 
converter  camera  shutter  operation  is  typical  of  system  operation  and 
will  be  discussed  in  some  detail.  Basic  time  information  was  derived 
from  the  Fastax  system.  Timing  marks  (placed  on  the  Fastax  film 
through  the  operation  of  a  light- chopper  driven  by  a  high  speed, 
precisely  synchronized  electric  motor)  permitted  accurate  determina¬ 
tion  of  Fastax  film  speed.  The  Fastax  film  distances  between  winkers 
marking  PM  tube  output  (Tq)  and  the  winkers  actuated  by  each  camera 
shutter  signal  were  then  converted  to  times.  On  the  Dynafax  film, 
distances  were  determined  between  the  winkers  and  the  winkers 
marking  camera  shutter  operation.  (Two  winkers  were  utilized,  one 
on  each  side  of  the  film  track  to  correct  for  possible  slit  misalign¬ 
ment).  Division  of  each  distance  on  the  Dynafax  film  by  its  appropriate 
time  from  the  Fastax  calculations  furnished  three  values  for  Dynafax 
film  speed.  The  average  of  the  three  values  was  then  adopted  as  the 

♦  Beckman  &  Whitley  "Dynafax"  Model  3193  camera. 
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true  Dynafax  film  speed.  The  distances  on  the  Dynafax  film  between 
the  impact  flash  and  the  camera  shutter- signal  winkers  were  then 
divided  by  Dynafax  film  speed  to  yield  the  times  from  impact  to  each 
camera  shutter  operation.  Raw  data  and  computations  for  shot  #2464 
appear  in  Appendix  B.  A  photograph  of  the  Dynafax  camera  and 
associated  winker  system  appears  as  Fig.  9. 

Accurate  determination  of  effective  interframe  times  in  the 
B&tW  300  camera  is  extremely  difficult.  Between  any  two  frames, 
errors  of  as  much  as  10%  of  the  nominal  (or  average)  interframe  time 
seem  probable.  Over  an  interval  of  many  frames,  the  resultant  error 
would  proba.bly  be  less  because  errors  would  tend  to  cancel  out.  Since 
the  exposure  time  is  roughly  25%  of  the  interframe  time,  the  precision 


Fig.  9.  Photograph  of  Dynafax  Camera 
and  Associated  Winker  System 
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of  the  data  acquired  is  limited  to  about  .2  *  j  .4  mm  in  object  plane 
spatial  resolution  and  .25  to  ,  50  microseconds  in  time  resolution, 

Switches 

Throughout  the  course  of  the  experiments,  various  types  of 
switches  were  employed  to  initiate  the  timing  cycle.  The  two  types 
employed  with  the  B&W  300  camera  were  mylar- foil  switches,  closed 
by  projectile  perforation  (Ref  2;93)  and  ion-probe  switches,  closed  by 
projectile  passage  between  the  probes  (Ref  28:18).  Fig.  10  presents  a 
schematic  signal-flow  diagram  of  the  instrumentation  discussed  thus 
far.  The  output  signal  of  the  initiating  switch  (ion-probe  or  mylar- 
foil)  was  directed  to  actuate  the  light  source,  either  directly  or  through 
a  time  delay  generator  depending  upon  the  distance  of  the  switch  from 
the  target  face.  The  initiating  signal  simultaneously  actuated  three 
separate  time  delay  generators.  The  output  of  each  generator  then 
fired  appropriate  x-rays.  The  last  time  delay  generator  simul¬ 
taneously  activated  the  Kerr  cell,  thereby  shutting  off  the  source  of 
light  for  the  B&W  300  camera.  The  initiation  or  T^  signal,  the  output 
of  the  three  x-ray  time  delay  generators,  and  the  firing  signal  from 
each  x-ray  head  switch  wer?  recorded  on  the  Fastax  film  rec  ord 
through  the  activation  of  xenon  winkers.  The  T^  signal  and  the  outputs 
of  the  first  and  last  time  delay  generators  were  likewise  recorded  on 
the  Dynafax  fJiin  record,  along  with  the  Impact  flash. 
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The  mylar-foil  and  ion-probe  initiating  switches  employed  in 
conjunction  with  the  B&tW  300  camera  did  not  prove  entirely  satisfac¬ 
tory  for  the  purposes  of  this  experiment.  The  mylar-foil  switch 
Invariably  functions  correctly:  but  once  closed,  It  frequently  gives 
repeated  output  signals  as  the  triggering  networks  repeatedly  discharge 
through  the  short-circuited  switch.  Also,  debris  resulting  from 
penetration  of  the  switch  proceeds  downrange  with  the  projectile  and 
interferes  with  the  desired  optical  record  of  the  ejecta  cloud  forma¬ 
tion.  While  the  ion-probe  switch  la  free  from  the  problems  of  debris 
and  repeated  output  signals,  it  is  a  much  less  positive  device.  For 
these  reasons,  it  was  decided  to  employ  a  photomultiplier  (PM)  tube 
to  sense  the  impact  flash.  Because  the  output  of  the  PM  tube  is  only 
a  few  volts,  it  was  necessary  to  amplify  the  output  to  a  level  sufficient 
to  activate  the  time  delay  and  thyratron  networks.  A  signal  flow 
diagram  of  the  resulting  system  appears  in  Fig.  11.  The  output  of 
the  signal  amplifier  was  recorded  on  the  Fastax  and  Dynafax  film 
records  along  with  the  impact  flash  (Dynafax  only)  and  image  converter 
camera  shutter  operations.  Recording  of  the  impact  flash  was 
necessary  since  measurements  indicated  that  there  was  a  delay  of 
about  1  microsecond  between  PM  tube  activation  and  the  resulting 
winker  firing. 

Target  Tank  Pressure 

Because  aerodynamic  drag  is  required  to  separate  from  the 
projectile  the  plastic  sabots  currently  in  use  at  the  AFML  facility,  it 
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was  necessary  to  maintain  range  pressure  no  lower  than  25.  2  torr  to 
ensure  satisfactory  separation.  Early  cloud  measurements  indicated 
that  even  this  slight  an  atmosphere  was  sufficient  to  decelerate 
appreciably  the  extremely  small  particles  constituting  the  early  ejecta. 
Therefore  for  the  image  converter  camera  experiments,  the  target 
tank  was  isolated  from  the  remainder  of  the  range  by  a  6.  5  ^  mylar 
diaphragm.  The  pressure  in  the  isolated  target  tank  was  then  reduced 
to  a  range  of  .  09-.  175  torr  by  means  of  a  supplementary  vacuum  pump. 
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Fig.  12.  Photograph  of  Microviewer  Used  for 
Analyzing  Optical  and  X-ray  Films 

Actual  crater  depth  and  diameter  measurements  were 
obtained  using  a  depth  gage  in  conjunction  with  the  calibrated-travel 
table  of  the  microviewer  (Fig.  13).  The  accuracy  of  measurement  of 
both  depth  and  diameter  was  on  the  order  of  t.  05  mm.  However, 
because  the  craters  are  not  perfectly  symmetrical,  the  values  deter¬ 
mined  for  crater  diameter  depend  to  some  extent  on  the  axis  of  meas¬ 
urement.  A  review  of  the  variation  in  measurements  (Appendix  B) 
indicates  that  the  determined  values  are  probably  accurate  to  about 
iO,  3  mm.  The  method  of  measurement  is  shown  in  Fig.  14.  Maximum 
depth  of  the  crater  (Pqjt)  was  measured  from  the  undisturbed  target 
face.  Crater  diameter  (D^-p)  was  established  as  the  mean  diameter 
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represented  by  the  intersection  of  the  elevation  of  the  undisturbed 
target  face  with  the  wall  of  the  crater. 


Projectiles  and  Tarpeta 

All  experiments  were  conducted  using  3.  ?.  mm  spheres  of  2017 
aluminum  as  projectiles.  All  targets  were  cut  from  a  single  ingot  of 
1100-0  aluminum.  The  ingot  was  specially  cast  and  an  attempt  made 
to  provide  grains  of  uniform  size.  Subsequent  metallurgical  and 
tensile  tests  showed  that  all  significant  properties  were  homogeneous 
with  respect  to  both  position  and  orientation.  The  ingot  possessed  a 
uniform  Brinell  hardness  value  of  24  (Ref  25:32),  Target  dimensions 
(in  cm)  were  as  follows: 


Shot  No. 

Length 

Height 

Depth 

2384 

8.9 

3.  6 

2.4 

2385 

00 

4.  3 

2.4 

2386 

8.  9 

5.0 

2.4 

2387 

7.  6 

5.  1 

3.  1 

2463 

8.  3 

7,  7 

4.  0 

2464 

8.  3 

7.  7 

4.2 

lU,  Crate?  Growth  Investigation 


X-ray  Layout 

In  any  hypervelocity  impact  on  an  opaque  target,  the  ejecta 
cloud  precludes  direct  optical  viewing  of  the  crater  formation  process. 
X-rays  capable  of  penetrating  the  target  material  provide  the  most 
accurate  means  of  following  crater  growth  in  these  situations.  A 
battery  of  five  sequential  flash  x-rays  were  used  in  this  investigation 
to  generate  data  on  crater  growth  as  a  function  of  time. 

In  order  to  obtain  profiles  of  both  crater  diameter  and  depth, 
two  300  kv  x-ray  heads  were  positioned  vertically  above  the  tax  get 
tank.  They  were  aimed  at  the  target  through  a  plexiglass  window  in 
the  top  of  the  tank.  The  x-ray  film  (Kodak  Royal  Blue  Medical)  was 
placed  in  a  cassette  directly  under  the  target  holder  in  the  bottom  of 
the  target  tank.  Experimental  Dupont  intensifier  screens  were  used 
in  the  cassette  to  enhance  the  contrast  obtainable  in  x-raying  through 
the  solid  aluminum  target.  A  schematic  of  the  vertical  x-ray  layout 
appears  in  Fig.  15. 

The  remaining  three  x-ray  heads  (one  300  kv  and  two  105  kv) 
were  positioned  at  the  rear  of  the  target  tank.  They  were  similarly 
aimed  at  the  target  through  a  plexiglass  window  in  the  rear  of  the  tank. 
The  x-ray  film  was  placed  in  a  cassette  mounted  vertically  in  the 
forward  portion  of  the  target  tank.  A  hole  was  bored  in  the  middle  of 
the  cassette  to  provide  a  means  of  projectile  passage,  Dupont 
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Fig.  15.  Positioning  of  Vertical  X-ray  Heads 
(Dimensions  are  In  cm) 


intensifier  screens  were  also  employed  in  this  cassette  to  enhance 
contract.  A  schematic  of  the  layov:..  .appears  in  Fig.  16. 

Since  x-rays  can  be  considered  neither  a  point  nor  a  parallel 
source,  determinations  of  tj’u '  •’iirensions  from  a  single  radiograph 
are  very  difficult.  However,  haiacteristic  dimensions  of  a 
symmetrically  changing  body  (such  as  a  crater)  can  be  readily  deter¬ 
mined  by  using  the  technique  of  comparing  the  measurements  on  two 
radiographs,  one  taken  during  the  event  and  the  other  after  the 
process  is  complete.  Of  course,  the  geometrical  relationships 
between  x-ray  head,  target,  and  film  plane  must  be  identical  for  the 
two  radiographs.  The  diameter  and  depth  of  the  crater  during  the 
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Fig.  16.  PoBitioning  of  Horizontal  X-ray  Heads 
(DimensionB  are  in  cm) 


event  are  then  readily  obtained  by  applying  the  ratio  of  measured  i 

dimensiona  on  the  two  radiographs  to  the  final  crater  parameter  as  j 

determined  by  actual  measurement  of  the  cratered  target.  1 

1 

1 

In  order  to  ensure  that  the  geometrical  relationships  could  be  j 

held  constant,  a  combined  target-  and  x-ray  cassette-holder  was  con-  ] 

i 

I 

structed  (Fig.  17).  The  target  was  supported  on  a  plexiglass  ledge  pro-  i 

1 

vided  with  a  lip  to  hold  the  target  securely  against  the  vertical  fiber-  j 

board  target-holder.  The  plexiglass  ledge  (with  target  resting  on  it) 
was  supported  on  bolts  extending  from  the  fiherboard  target-holder. 

Wing  nuts  were  affixed  to  the  bolts  to  facilitate  target  mounting  and 
removal  and  to  allow  the  target  to  be  held  securely  against  the  holder. 

I 

I 
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To  prevent  multipie-exposure  of  the  x-ray  film  resulting  from 
the  firing  of  the  three  separate  horirontal  heads,  a  6.4  mm  lead 
shield  was  placed  behind  and  bolted  to  the  target-holder.  A  4,  76  cm 
hole  was  then  bored  through  both  lead  shield  and  fiberboard  target- 
holder  directly  behind  the  target  position.  Thus  maximum  x-ray 
energy  from  each  head  was  allowed  to  pass  through  the  target  and 
strike  different  locations  on  the  film  in  the  vertical  cassette. 

In  order  to  ensure  that  the  x-ray  cassette  could  be  relocated 
for  the  radiograph  to  be  taken  after  the  shot  in  the  same  position  it 
had  occupied  for  the  radiograph  taken  during  crater  growth,  an  x-ray 
cassette  holder  was  similarly  constructed  of  fiberboard.  Plexiglass 
ledges  were  bolted  to  the  fiberboard  at  the  top  and  bottom,  spaced  so 
that  the  cassette  fit  snugly  between.  Thumb  screws  were  added  to 
clamp  the  cassette  in  position  on  the  ledges.  Through  the  2.  54  cm 
hole  bored  in  both  cassette  o.nd  fiberboard  was  passed  a  short  length  of 
pipe.  This  assured  correct  horizontal  positioning  and  helped  to  make 
the  x-ray  cassette  light-tight  around  the  hole.  (The  area  where  the 
pipe  emerged  from  the  cassette  was  repacked  with  Duxseal*  each  tirne 
the  cassette  was  loaded  with  film.  )  Finally,  both  target-holder  and 
x-ray  cassette  holder  were  rigidly  fixed  at  the  desired  spacing  by 
parallel  sections  of  1.27  cm  diameter  threaded  rod  at  the  top  and 
bottom..  The  entire  unit  was  affixed  to  the  target  tank  by  the  same 

*  Brand  name  for  a  putty-like  sealing  and  caulking  compound  used  In 
the  electrical  iridustry. 


36 


threaded  rod.  A  photograph  of  the  completed  unit  appears  In  Fig.  18. 

The  horizontal  x-ray  cassette  (containing  film  to  be  exposed  by  the  two  1 

vertically-mounted  300  kv  heads)  was  supported  in  the  bottom  of  the 
target  tank  on  lead  bricks  and  w'edged  securely  under  the  target 
holder. 

Radiographic  Re  suits 

The  radiographs  from  shot  #2386  appear  as  Figs.  19  and  20. 

(Radiographs  of  shots  #2384,  2383,  and  2387  appear  as  Figs.  37 
through  42  in  Appendix  B.  )  To  the  naked  eye,  the  contrast  is  generally 
adequate  to  determine  crater  limits.  However,  the  magnification 
resulting  from  viewing  the  negatives  projected  through  the  lOX  head 


Fig.  18.  Photograph  of  Target  and  Vertical 
X-ray  Cassette  Holders 
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During  Shot 


After  Shot 


13.48^  sec 


15.66^ 


stc 


Fig.  19.  Radiographs  of  Crater  Profile  During  and  After 
Shot  #2386  (Vertical  X-ray  Setup) 
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During  Shot 


After  Shot 


14. 49  u  sec 


Fig.  20.  Radiographs  of  Crater  Diameter  During  and  After 
Shot  #2386  (Horizontal  X-ray  Setup) 
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of  the  microviewer  greatly  diminished  the  contrast.  In  order  to 
provide  a  guide  to  the  location  of  the  crater  limits,  light  pin-pricks 


were  made  along  the  crater  limits.  When  placed  on  the  microviewer, 
the  pin-pricks  made  it  possible  to  determine  crater  limits.  On  the 
profile-view  pictures,  crater  diameter  was  measured  along  the  line  of 
the  undistorted  target  face;  maxim.um  depth  was  similarly  measured 
from  this  line.  The  angle  between  a  tangent  to  the  crater  wall  at  the 
point  or  original  target  surface  and  the  target  face  was  determined  for 
each  profile  view.  The  angle  remained  virtually  constant  (Fig.  21). 

I 

On  the  radiographs  taken  from  behind  the  target,  the  maximum  width  I 

across  the  crater  image  on  a  horizontal  line  was  measued.  Identical 
procedures  were  used  on  both  the  radiographs  taken  during  and  those 
taken  after  the  shot. 

Film  measurements  of  crater  depth  and  diameter  appear  in 
Table  II.  From  the  film  measurements,  the  ratios  of  crater  para¬ 
meters  during  and  after  the  shot  were  computed.  The  resultant  time- 
history  of  crater  growth  is  presented  graphically  in  Figs.  22  and  23. 

Since  the  growth  of  both  crater  diameter  and  depth  as  deter¬ 
mined  in  this  study  appeared  to  follow  an  exponential  growth  law,  a 
curve  of  the  form  1-e"^^,  where  t  is  measured  in  microseconds,  was 
fit  to  the  data.  Best-fits  were  obtained  for  values  of  k  =  .220  in  the 
case  of  diameter  growth  (standard  deviation  =  4.  9%  of  final  crater  dia¬ 
meter)  and  k  =  .  2  14  in  the  case  of  depth  growth  (standard  deviation  = 

2.83%  of  final  crater  depth).  In  the  latter  case,  using  the  value  of 
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Kig.  21.  Angle  Between  Tangent  to  Crater  Wall  at  Original  Target  Surface  and  Target  Face 


ray  Film  Measurement*  of  Crater  Diameter  and  Depth 
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Crater  Diameter  Growth 
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Fig.  23.  Plot  of  Crater  Depth  Growth 


k  =  .  220  as  determined  for  diameter  growth  provides  a  surprisingly 
good  fit,  with  a  standard  deviation  of  only  2.  96%.  The  best-fit  curves 
are  plotted  on  Figs.  22  and  23. 

Measurements  of  Final  Craters 

A  photograph  of  the  six  cratered  targets  appears  as  Fig.  24. 
Since  the  velocities  of  impact  and  projectile  weights  were  controlled 
as  carefully  as  possible  to  ensure  constant  impact  energy,  the  result¬ 
ing  craters  are  very  similar,  as  might  be  expected.  Spallation,  clearly 
visible  on  the  target  from  shot  #2384,  did  not  result  in  any  significant 
alteration  of  crater  size  or  shape.  Maximum  crater  depth,  measured 
from  the  plane  of  the  undisturbed  target  face,  varied  by  only  .  66  mm. 
Two  separate  diameter  measurements  along  perpendicular  traverses 
were  made  for  each  crater,  the  diameter  being  measured  from  the 
intersection  of  the  crater  wall  with  the  plane  of  the  undisturbed  target 
face  in  each  cas,e.  The  maximum  variation  in  crater  diameter  was 
,  7  mm.  The  resulting  values  of  crater  diameter  and  depth  are  pre¬ 
sented  in  Table  III.  (Raw  data  and  method  of  reduction  appear  in 
Appendix  B.  ) 

In  an  effort  to  determine  the  degree  of  reproducibility  of  the 
experiments,  the  crater  diameters  were  scaled  to  shot  #2386  as  a 
standard  according  to  the  cube  roots  of  impact  energy.  /Both  Riney 
and  Walsh  concur  in  energy  scaling,  the  latter  with  a  slight  reservation 
(Ref  3;257-258Jy^  As  can  be  seen  from  Tabic  IV,  the  resulting 


Shot  #2384 


Shot  #238S 


Shot  #2386 


Shot  #2387  Shot  #2463  Shot  #2464 

Fip,  24.  Photographs  of  Cratered  Targets 


4 


TABLE  III 

Direct  Measurements  of  Crater  Diameter  and  Depth 


Shot 

Numbe  r 

Average  Crater 
Diameter 
(cm) 

Maximum  C 
Depth 
(cm) 

2384 

1. 65 

.  90 

2385 

1, 67 

.  93 

2386 

1. 64 

.  92 

2337 

1.  60 

CO 

2463 

i  .  03 

.  93 

2464 

1.  60 

.  87 

TABLE  IV 
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variations  from  standard  are  of  the  same  order  as  the  variations  in 


uncorrected  crater  diameters.  Hence,  it  is  felt  that  the  variations 
represent  normal  measurement  scatter,  and  that  the  craters  are 
virtually  identical.  It  would  therefore  appear  that  the  experiment  is 
entirely  reproducible. 
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Fig.  32.  Examples  of  Late-Time  Framing  Camera  Photograph 
of  Ejecta  Cloud  From  Shot  f2387  With  Times  From. 
Impact  sec)  ^Magnification  =  0.  14  (non'jy 


™  '  -{A 


IV.  EjectA  dead  lnv«*8tigation 


Camera  Layouts 

Franiing  Camera.  Because  of  the  rapidity  with  which  tlie  ejecta 
cloud  is  formed,  ultrahigh  speed  optical  photography  provides  the  best 
tool  for  observing  its  formation  and  growth,  A  Beckman  &  Whitley 
Model  300  framing  camera  v/as  initially'  used  to  obtain  optical  photo¬ 
graphs  of  the  ejecta  cloud  (see  Fig.  25).  At  the  extremely  high  framing 
rates  (4.5  million  frames/sec)  which  this  camera  is  capable  of, 
shuttering  is  best  accomplished  by  accurately  controlling  the  synchron¬ 
ization  and  duration  of  the  back  lighting  source.  Since  the  spark 
light  source  em.ployed  had  an  extremely  rapid  rise  time  (on  the  order 
of  one  microsecond),  it  was  necessary  only  to  control  light  extinction. 
This  was  accomplished  by  using  a  Kerr  cell  which  effectively  cut  the 
light  off  in  less  than  0.  1  microseconds.  A  schematic  of  tlit  >tical 
layout  employed  with  the  B&W  300  caniera  appears  in  Fig.  26.  An 
integrated  schematic  of  both  optical  and  x-ray  layouts  appears  in 
Fig.  27. 

The  large  ratio  of  frame  exposure  time  to  interframe  time 
(1:4),  low  optical  resolution,  and  inaccuracies  in  interframe  time  com¬ 
bined  to  reduce  the  precision  of  the  resultant  data  to  values  unaccept¬ 
able  for  the  cloud-edge  investigation.  However,  data  from  this  camera 
proved  sufficient  for  monitoring  gross  cloud  dimensions  and  the  motion 


of  discrete  ejecta  particles. 


Camera  Layouts 


Ejecta  Cl  Pad  Investigation 


Framing  Camera.  Because  of  the  rapidity  with  which  the  ejecta 
cloud  is  formed,  ultrahigh  speed  optical  photography  provides  the  best 
tool  for  observing  its  formation  and  growth.  A  Beckman  &  Whitley 
Model  300  framing  camera  v/as  initially  used  to  obtain  optical  photo¬ 
graphs  of  the  ejecta  cloud  (see  Fig.  25).  At  the  extremely  high  framing 
rates  (4.  b  million  frames/sec)  which  this  camera  is  capable  of, 
shuttering  is  best  accomplished  by  accurately  controlling  the  synchron¬ 
ization  and  duration  of  the  back  lighting  source.  Since  the  spark 
light  source  employed  had  an  extremely  rapid  rise  time  (on  the  order 
of  one  microsecond),  it  was  necessary  only  to  control  light  extinction. 
This  was  accomplished  by  using  a  Kerr  cell  which  effectively  cut  the 
light  off  in  less  than  0.  1  microseconds.  A  schematic  of  tht  >tical 
layout  employed  with  the  B&W  3C0  camera  appears  in  Fig.  26,  An 
integrated  schematic  of  both,  optical  and  x-ray  layouts  appears  in 
Fig.  27. 

The  large  ratio  of  frame  exposure  time  to  interframe  time 
(1;4),  low  optical  resolution,  and  inaccuracies  in  interframe  time  com¬ 
bined  to  reduce  the  precision  of  the  resultant  data  to  values  unaccept¬ 
able  for  the  cloud-edge  investigation.  However,  data  from  this  camera 
proved  sufficient  for  monitoring  gross  cloud  dimensions  and  the  motion 
of  discrete  ejecta  particles. 
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Fig.  25.  Photograph  of  Beckman  &  Wliitley  Model  300  Framing 
Camera  With  Example  of  Picture  Format 


LENS 


Fig.  26.  Schematic  of  Arrangement  of  Beckman  &  Whitley  Model  300 
Framing  Camera  With  Light  Source  and  Kerr  Cell 


VERT  I 


Integrated  Schematic  of  Instrumentation  System  Used  in 
Experiments  Employing  Framing  Camera  and  X-ray  Setups 


Image  Converter  Camera  System.  A  three-frame  image 


converter  camera  system  was  constructed  to  replace  the  framing 
camera  as  the  principal  source  of  cloud- edge  data.  The  extremely 
short  exposure  times  of  the  individual  cameras  (10  ns*  and  1  ns**) 
together  with  variable  and  precisely  measured  interframe  times  pro¬ 
vided  an  accurate  time  base  for  the  cloud- edge  experiments.  The  high 
and  uniform  optical  resolution  of  the  image  converter  cameras  (15  line 
pairs  per  mm  on  the  film  plane)  was  sufficient  to  resolve  all  elements 
of  the  ejecta  cloud  of  interest  to  this  investigation. 

The  image  converter  cameras  were  aligned  on  the  same  optical 
arts  thTO’igh  a  system  of  mirrors  as  shown  in  Fig.  28.  The  topmost 
(1  ns)  camera  viewed  the  target  through  an  angled  mirror,  silvered  to 
transmit  10%  and  reflect  90%  of  the  incident  light.  The  second  or 
middle  camera  was  positioned  below  the  first.  It  viewed  the  target  by 
reflection  off  both  the  top  mirror  and  a  second  mirror,  parallel  to  and 
directly  below  the  first,  silvered  to  transmit  60%  and  reflect  40%. 

The  bottom  camera  was  aimed  vertically  upward  and  positioned 
directly  beneath  both  mirrors.  Hence  it  viewed  the  target  through 
the  lower  mirror  and  by  reflection  off  the  upper  mirror.  A  photograph 
of  the  cameras  and  mirrors  in  position  appears  as  Fig.  29. 

*  Beckman  &  Whitley  Model  500  Image  Converter  Camera. 

**  Experimental  image  converter  camera  loaned  to  AFML  by 
Beckman  &  Whitley  Corporation. 
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Fif;.  29.  Photoj’raph  of  linage  Converter  Cameras 
and  Mirrors  in  Position  on  the  Range 
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Optical  Results 


Framing  Camera  Photographs.  From  photographs  of  a  grid,  it 
was  determined  that  the  B&W  framing  camera  introduced  a  slight 
amount  of  optical  distortion.  A  computer  program  was  written  to 
remove  this  distortion.  In  addition,  the  program  corrected  film 
dimensions  to  real  space  dimensions  by  application  of  a  magnification 
factor.  The  magnification  factor  was  determined  independently  for 
each  frame.  The  film  was  aligned  in  the  film  reader  so  that  the  y-axis 
was  parallel  to  the  vertical  arm  of  a  set  of  crosshairs  mounted  per¬ 
manently  in  the  camera.  Then  the  x-  and  y-coordinates  of  the  cross¬ 
hair  intersection  and  of  two  points  at  known  distances  apart  (in  real 
space)  were  read.  The  program  then  computed  the  appropriate  magni¬ 
fication  factor  to  reduce  all  x-y  film  positions  to  real  space  coordinates. 
A  simple  orthogonal  transformation  was  then  made  to  shift  the 
coordinate  origin  to  the  desired  center  of  impact  on  the  target  face. 

The  ejecta  cloud  photographs  of  shot  #2385  covered  a  time 
span  from  approximately  .79  to  3.45  ^  sec  after  impact.  Representa¬ 
tive  samples  of  these  photographs  appear  as  Fig.  30.  Upon  close 
observations,  it  can  be  seen  that  the  cloud  is  not  perfectly  symmetri¬ 
cal.  Hence  the  two  sides  were  treated  separately.  The  coordinates  of 
points  along  the  edge  of  each  side  of  the  cloud  in  each  frame  were  read 
on  the  microviewer  and  reduced  to  real- space  coordinates  by  the 
computer  program.  (A  length  of  plastic  threaded  rod  had  been  screwed 
into  the  target  holder  ledge  directly  under  the  target;  the  thread  ends 
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provided  the  required  kiiown-distance  measurement,  )  Polynomials 
were  then  fit  by  the  method  of  least  squares  to  the  cloud- edge 
coordinates.  A  satisfactory  fit  was  achieved  in  all  instances.  The 
maximum  RMS  error  for  the  top  cloud  edge  curve  was  1.47  mm;  for 
the  bottom  curve,  1.83  mm. 

Despite  this  apparent  satisfactory  fit  of  polynomials  to  thb 
cloud  edges,  the  cloud-edge  program  was  able  to  find  the  required 
straight  line  trajectories  for  only  a  very  limited  number  of  points  on 
these  curves.  A  summary  of  ‘he  computer  program  output  is  pre¬ 
sented  in  Table  V. 

It  is  quite  apparent  from  these  data  that  the  program,  with 
input  data  from  shot  #2385,  is  not  giving  realistic  outputs.  While  the 
origin  (and  presumably  the  crater  radius)  is  seen  generally  to  increase 
with  time,  the  entry  for  230  nanoseconds  reverses  this  trend.  The 
apparent  decrease  in  ejection  angle  with  time,  while  unexpected,  is 
not  impossible  in  light  of  Gault’s  previously  cited  findings.  However, 
for  the  ejection  velocity  to  increase  monoto;nically  with  time  is  a 
virUial  impossibility.  Were  this  the  case,  the  entire  model  would  be 
invalidated,  since  interparticle  collisions  would  almost  inevitably 
occur,  and  the  assumption  of  constant  velocity  after  ejection  would  no 
longer  be  tenable.  Lastly,  the  inconsistency  of  the  results  obtained  by 
changing  only  the  third  curve  (from  frame  22  to  frame  20)  strongly 
suggests  that  the  data  are  not  meaningful. 


TABLE  V 


Cloud-£dge 

Computer 

Program  Outputs, 

Shot  #2385 

Time  of 
Ejection 
(ns) 

Origin 

(cm) 

Angle 

(Deg) 

V  elocity 
(km/ sec) 

20 

.  285 

63.  4 

3.  61 

50 

lZI 

62.  1 

3.  87 

80 

.  288 

60.  6 

4.  17 

110 

.291 

59.  0 

4.  50 

140 

.294 

57.  3 

4.  88 

170 

.297 

55.  1 

5.  40 

220* 

,  253 

64.  5 

2.  1 1 

230 

.289 

51.  3 

6.  36 

350* 

.260 

65.  5 

2.  28 

*  Data  obtained  irom  frames  12,  16,  and  20.  All 
other  data  obtained  from  frames  12,  '6,  and  22. 


A  comparison  of  the  artificial  constant-velocity  curve  at  a  time 
after  impact  of  3.  58^  sec  ■with  the  actual  ejecta- edge  curve  of 
shot  #2335  at  3.45^  sec  is  presented  in  Fig.  31.  The  t'wo  cloud  edges 
appear  to  be  at  least  superficially  similar  in  their  lower  portions. 
However,  in  the  upper  portions,  the  actual  ejecta  cloud  appears  to  be 
flattening  out  more  rapidly  than  the  constant-velocity  cloud.  This 
suggests  the  possibility  that  the  early  ejecta  is  being  decelerated  by 
atmospheric  drag.  The  earliest  portions  of  the  ejecta,  being  the  least 
massive  and  the  fastest  would  naturally  undergo  the  greatest  decelera¬ 
tion. 


Fig,  31.  Comparison  of  Observed  Cloud  Edge  With 
Computer-Produced  Constant- Velocity  and 
Decelerating  Particle  Curves 

c  o 

■J  y 


To  teat  further  this  possibility,  a  deceleration  term  was  intro¬ 
duced  into  the  cloud  formation  computer  program.  Particle  velocity 
was  reduced  by  a  factor  of  e“^^,  where  t  =  the  time  of  flight  of  the 
particular  particle,  and  k  was  arbitrailly  selected  as  3  x  10“^  divided 
by  sequential  particle  number.  This  term,  then,  was  based  on  the 
assumption  that  the  earlier  particles  would  experience  more  decelera¬ 
tion  than  the  later,  more  massive  particles.  The  curve  resulting  from 
this  exercise  is  plotted  on  Fig.  31  as  a  dashed  line.  It  is  indistin¬ 
guishable  from  the  constant  velocity  computer  curve  in  its  lower 
region,  but  flattens  out  in  itc  upper  portion  In  generally  the  same 
fashion  as  does  the  actual  curve. 

The  framing  camera  photographs  from  shot  #2387  covered  the 
entire  time  span  from  impact  »\ntil  ejection  had  virtually  ceased. 
However,  because  the  Kerr  cell  did  not  operate  properly  to  cut  off  the 
light  source,  the  frames  are  rewritten,  1.  e. ,  in  about  half  the  frames, 
two  separate  images,  one  from  early  time  and  one  from  late  time,  are 
superimposed  in  the  same  frame.  The  early-time  frames  were 
deemed  unusable,  since  cloud  edges  could  not  be  accurately  deter¬ 
mined.  On  the  late-ti:r. .  frames  (Fig.  32),  discrete  particles  are 
readily  distinguishable,  A  number  of  these  particles  were  identified 
by  characteristic  shapes,  and  their  x-y  coordinates  determined  in  as 
many  successive  frames  as  was  possible.  The  film  positions  were 
then  converted  to  real-space  positions  through  the  computer  program 
previously  discussed.  By  dividing  the  x-  and  y-distance  traveled  by 
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I  Fig-  32.  Examples  of  Late-Tirp.e  Franiin 

I  of  Ejecta  Cloud  From  Shot  #'238' 

Impact  (u  sec)  ^^Magnification  = 


the  time  difference  between  frames,  x-  and  y-components  of  velocity 
(V^  and  Vy)  were  determined  for  each  particle.  Actual  particle  velocity 
(speed  and  axigle)  were  then  determined.  Dividing  the  final  x-distance 
of  the  particle  by  yielded  its  time  of  flight,  based  on  the  assump¬ 
tion  that  its  constant  velocity  free-flight  phase  had  originated  at  the 
target  surface.  Similarly,  this  time-of-flight  multiplied  by  gave 
its  y-travel,  and  this  distance  subtracted  from  its  final  y-distance 
yielded  the  assumed  radial  distance  of  particle  departure  from  the 
target  face. 

The  results  of  the  application  of  these  techniques  to  20  particles 
identifiable  in  frames  36  and  44  of  shot  #2387  are  presented  in 
Table  VI,  (Raw  data  and  sample  calculations  appear  in  Appendix  B.) 

The  data  are  listed  by  apparent  time  of  ejection  of  the  particle.  The 
negative  times  are  of  course  physically  impossible.  Since  the  time 
between  frames  of  the  B&W  300  camera  is  extremely  difficult  to  deter¬ 
mine  accurately,  it  is  entirely  possible  that  all  times  shown  could  be 
in  error  by  as  much  as  .  5  microseconds.  Whatever  time  error  exists 
would  be  applied  equally  to  all  particles. 

These  data  are  plotted  graphically  in  Fig.  33.  Clearly  the 
scatter  of  the  data  for  angle  and  velocity  preclv^es  determination  of 
any  trend.  In  the  case  of  the  origin  (or  assumed  crater  radius),  if 
the  data  points  for  -1.  31  and  -4.  38  microseconds  are  discarded,  the 
remaining  points  do  suggest  increasing  crater  radius  with  time. 
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TABLE  VI 


Crater  Parameters  Determined  from  Discrete 
Ejecta  Particle  Movement, 

Frames  36  to  44,  Shot  fi^2387 


Time  of 
Ejection 
(/I  Bee) 

Origin 

(cm) 

Angle 

(Deg) 

Velocity 
(km/ sec) 

00 

1 

.  551 

53.  8 

.  66 

-1.  31 

.  739 

57.  1 

.  67 

-1, 28 

.203 

57.  1 

1. 07 

1 

o 

.  320 

57.  5 

1.26 

1 

o 

o 

00 

.  381 

43.  8 

.  58 

0.  21 

.  386 

57,  4 

1. 08 

0.  50 

,406 

49.  9 

.  64 

0.  87 

.488 

59.  5 

.  98 

0.  93 

.  325 

55.  2 

.  94 

1.  30 

.  310 

54.  6 

.  91 

1. 42 

.  399 

44.  4 

.  57 

1 .  50 

,424 

61. 2 

1.  34 

1.  58 

.  381 

56,  2 

.  97 

1.  93 

.  381 

56.  0 

.  98 

1.  99 

.457 

60.  2 

1.  16 

2. 06 

.439 

54.4 

.  90 

2.  30 

.485 

57.4 

.  94 

2,  43 

.559 

55.  8 

.  79 

2.84 

.597 

52.  0 

.  64 

2,  85 

.663 

61. 6 

1.  36 
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PARTICLE  ORIGIN  (cm)  ANGLE  OF  EJECTION  (0E6)  PARTICLE  VELOCITY  {fcm/»«c) 


Fig.  33.  Crater  Parameters  vs.  Time  Determined  From  Analysis 
of  Discrete  Particle  Motion,  Shot  #2387 
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However,  the  scatter  of  even  the  remaining  data  is  such  as  to  vi'arrant 
questioning  its  validity. 

Image  Converter  Photographs.  In  an  effort  to  overcome  the 
problem  of  determining  precisely  photograph  times  with  respect  both 
to  impact  and  to  each  other,  three  optically  aligned  Image  converter 
cameras  were  used  to  photograph  the  ejecta  clouds  of  shots  #2463  and 
#2464.  For  these  shots  target  tank  pressure  was  reduced  as  low  as 
possible  (.09-.  175  torr).  Spatial  oriiJiitation  was  achieved  by  photo¬ 
graphing  a  grid  aligned  with  the  projectile  trajectory  after  each  shot. 
Prior  to  target  removal,  the  exact  pomt  of  impact  was  determined  by 
sighting  down  the  gun  tube  with  a  telescope.  The  target  was  then 
removed  from  the  target  holder  and  the  grid  positioned  so  that  the 
intersection  of  two  specified  grid  lines  marked  the  point  of  projectile 
impact.  This  was  taken  as  the  origin  of  the  real- space  x-y  coordinate 
sy  stem. 

Measurements  of  the  grid  photographs  on  the  microviewer 
indicated  that  optical  distortion  was  negligible  (on  the  order  of  1/4%). 
Therefore  a  single  multiplication  factor  for  each  im.age  converter 
photograph  was  determined  from  the  grid  photographs.  The  x-  and 
y~axes  from  each  grid  photograph  were  transferred  to  the  appropriate 
actual  shot  photograph  by  scribe  marks,  made  while  viewing  the 
superimposed  negatives  on  a  light  table.  Examples  of  the  image 
converter  photographs  appear  in  Fig.  34. 
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An  example  of  raw  data  from  the  Image  converter  photographs 
of  shot  #2464  appears  in  Appendix  D,  As  in  the  case  of  the  BSiW  300 
ejecta  cloud-edge  photographs,  polynomials  were  fit  to  the  real-space 
coordinates  of  the  cloud  edges  determined  fron-.  the  image,  converter 
photographs  using  a  standard  least- squares  program.  The  two  edges 
of  the  cloud  were  treated  separately  because,  as  in  the  earlier  case, 
the  cloud  edges  v/ere  not  synimetrical.  Fourth-order  polynomials  were 
employed  in  all  cases,  since  the  maximum  RMS  error  of  0.  3  mm  was 
deemed  acceptable.  Examples  of  the  computer  progran^,  output  appear 
in  Table  VII. 

Unlike  the  earlier  results,  the  computer  output  re.sulting  from 
the  image  converter  photograph  data  appears  internally  consistent, 
i.  e.  ,  origin  of  ejection,  angle,  and  velocity  all  appear  to  be  behaving 
reasonably  within  the  range  of  about  4,42  to  5.69  microseconds  after 
im.pact.  Since  the  time  of  the  first  photograph  was  only  8.  52  micro¬ 
seconds  after  impact,  erratic  outputs  for  times  approaching  this  value 
are  to  be  expected  due  to  the  very  short  trajectory  distance  between 
the  "A"  curve  and  the  origin.  However,  the  failure  of  the  program  to 
obtain  data  for  times  earlier  than  4.42  microseconds  is  not  as  easily 
understood.  Analysis  of  the  reason  for  convergence  failure  within  the 
program  indicates  that  just  as  was  the  case  with  shot  #2385,  the  third 
or  "C"  curve  is  flattening  out  to  such  an  extent  that  a  correct  geo¬ 
metric  straight-line,  constant-velocity  solution  cannot  be  obtained. 

As  can  be  inferred  from  Fig.  31,  particle  deceleration  appears 
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TABLE  VU 


Cloud-Edge  Computer  Program  Outputs,  Shot  #2464 


Time  of 
Ejection 
aec) 

Origin 

(cm) 

Angle 

(Deg) 

Velocity 

(km/sec) 

4.  42 

.282 

61.  5 

3.  42 

4.  54 

.  309 

62.2 

3.42 

4.  67 

.  324 

62.4 

3.  41 

4.  79 

.  339 

62.  6 

3.  40 

4.  91 

.  353 

62.  7 

3.  39 

5.  02 

.  367 

62.  9 

3.  37 

5.  16 

.  384 

63.  1 

3.  35 

5.27 

.  397 

63.  2 

3.  33 

5.  37 

.410 

63.4 

3.  30 

5.  47 

.423 

63.  5 

3.  26 

5.  55 

.435 

63.  6 

3.  20 

5.  63 

.  446 

63.  6 

3.  14 

5.  68 

.  457 

63.  6 

3.  04 

5,  69 

.463 

63.  6 

2.  97 

t>8 


capable  of  producing  this  effect.  A  further  reduction  in  target  tank 
pressure  might  prevent  particle  deceleration,  but  such  vacuums  are 
difficult  to  secure  and  maintain. 

Although  neither  the  cloud-edge  computer  program  nor  the 

discrete-particle  computations  were  successful  in  consistently  deter- 

/' 

mining  crater  parameters  as  functions  of  time,  the  changes  in  width 
of  the  base  and  throat  of  the  ejecta  cloud  present  an  interesting  insight 
into  the  method  of  crater  formation.  These  two  cloud  parameters 
were  determined  in  as  many  photographs  as  possible,  and  the 
resulting  data  are  plotted  on  Fig.  35  along  with  the  best-fit  curve  of 
crater  diameter  growth.  The  data  are  presented  as  ratios  of  the  value 
of  the  appropriate  parameter  to  its  final  value,  as  determined  from 
late-time  photographs  in  which  the  parameters  stopped  increasing  and 
exhibited  only  random  variations  around  the  (assumed)  final  value.  In 
very  early  times,  the  base  width  increases  more  rapidly  than  the 
actual  crater  as  can  be  seen  from  Fig.  35.  The  throat  lags  behind 
both,  clearly  indicating  that,  in  these  times,  the  crater  is  expanding 
faster  than  the  radial  component  of  effective  ejecta  velocity. 


V .  Discussion  of  Results ,  Conclusions ,  and  Recommendations 


Discussion  of  Results 

The  purpose,  instrumentation,  and  general  results  of  each  of 
the  six  shots  contributing  data  to  this  investigation  are  summarized  in 
Table  VIII,  The  internal  consistency  of  crater  growth  data  (Figs.  22 
and  23)  acquired  from  the  radiographs  is  excellent.  The  shape  of  the 
diameter  growth  curve  (Fig.  22)  is  in  general  agreement  with  results 
obtained  by  Charters  (Ref  3;281)  for  1100-0  aluminum-on- 

aluminum  impacts  at  only  slightly  higher  velocities  (7.4  km/sec  vs. 

7.  0  km/sec  used  for  this  study).  However,  the  absolute  diameters 
reported  by  Charters  are  30-40%  greater  than  the  crater  diameters 
achieved  in  this  study.  /Crater  parameters  measured  during  this 
study  appear  consistent  with  previous  AFML  results  and  with  other 
published  results  (Ref  26)^  The  cause  of  the  disagreement  is  felt  to 
lie  in  Charters'  measurement  technique.  Apparently  he  measured  the 
base  of  the  ejecta  cloud  and  considered  it  equivalent  to  crater  diameter 
(Ref  3:280).  Results  achieved  in  the  study  being  reported  show  that 
there  is  an  appreciable  and  varying  difference  between  these  two 
parameters  (see  Fig.  22).  To  further  clarify  the  relationships  between 
cloud  and  crater  dimensions,  the  ejecta  cloud  base  diameters,  throat 
widths,  and  crater  diameter  growth  curve  of  the  present  study  are 
plotted  in  Fig.  36  together  with  Charters'  curve  of  crater  diameter 
growth.  It  is  doubtful  that  the  slight  difference  in  velocities  between 
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Table  VIII  (Continued) 
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Fig.  36.  Comparison  of  Ejecta  Cloud  Base  and  Throat  Data  With  Crater  Diameter  Growth  Curves 


the  two  experiments  could  account  for  the  large  disparity  in  reported 
diameters.  In  any  event,  the  figure  demonstrates  that  cloud  base 
cannot  be  considered  equivalent  to  actual  crater  diameter. 

Both  crater  diameter  and  depth  data  were  found  to  follow  very 
closely  a  growth  law  of  the  form  R  =  1-e”^^,  where  R  is  the  ratio  of 
transient  to  final  crater  dimension,  k  is  the  fitting  parameter 
(k  0.  22  for  the  case  considered),  and  t  =  time  in  microseconds  after 
impact.  The  decision  to  apply  such  an  exponential  curve  to  the  data 
was  based  solely  on  the  appearance  of  the  data,  and  no  physical  justi¬ 
fication  for  this  type  of  growth  law  is  proposed.  If  such  an  exponential 
growth  law  were  to  prove  applicable  over  wider  etiergy  ranges  and  for 
a  variety  of  target  materials,  it  would  provide  a  valuable  tool  for 
testing  theoretical  treatments. 

The  radiographs  demonstrate  that  the  crater  assumes  its 
characteristic  shape  very  early  in  the  cratering  process  and  retains 
this  shape  as  it  grows  in  size.  Ideal  hypervelocity  impact  craters  are 
almost  universally  described  as  "hemispherical".  This  judgment 
appears  to  be  based  largely  (if  not  entirely)  on  the  close  agreement 
between  values  of  depth  and  radius.  But  the  tangent  to  a  hemispherical 
crater  at  the  target  surface  would  always  make  an  angle  of  90  degrees 
with  the  target  face.  The  value  of  this  angle  determined  in  this  study 
remained  almost  constant  at  77  i 4  degrees.  It  is  therefore  concluded 
that  the  true  shape  of  the  crater  is  not  hemispherical  but  is  better 
represented  by  a  paraboloid. 
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The  x-ray  techniques  employed  in  this  study  appear  very 
promising.  In  particular,  the  comparative  technique  for  obtaining 
crater  measurements  appears  to  offer  advantages  of  ease  and  sim¬ 
plicity  over  any  absolute  method,  such  as  the  use  of  reference  grids. 
The  capability  of  making  usable  radiographs  c  transient  craters  in 
solid  aluminum  targets  has  been  demonstrated.  It  is  felt  that  careful 
analysis  of  radiographs  obtained  earlier  in  the  process,  when  the 
crater  lip  is  first  being  formed,  may  give  additional  insight  into  the 
mechanics  of  crater  formation  and  ejection. 

The  fact  that  the  crater  achieves  its  final  shape  early  in  the 
crater  growth  process  increases  the  problems  associated  with  using 
ejecta  cloud  information  to  predict  crater  growth  parameters.  It  was 
discovered  that  no  reasonable  estimate  of  crater  size  and  shape  could 
be  obtained  from  the  ejecta  cloud  photographs.  The  failure  of  the 
computer  program  to  predict  correctly  the  time  variation  of  crater 
grow'th  parameters  must  result  from  one  or  more  of  the  following 
causes:  (1)  inaccuracies  in  determining  true  ejecta  cloud  edge  position 
from  the  photographs,  (2)  inaccuracies  in  determ.ining  the  times  at 
which  the  cloud  edge  photographs  were  taken,  or  (3)  invalidity  of  the 
assumption  that  the  ejecta  travels  in  straight  lines  at  constant 
velocities. 

The  importance  of  timing  and  positional  errors  was  investi¬ 
gated  by  a  series  of  program  sensitivity  tests.  Artificially  generated 
cloud  curves  (described  in  Chapter  II)  were  altered  by  flattening  out 
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the  "C  curve.  The  resulting  curve  resembled  the  curve  produced  by 
decelerating  particles.  Then  both  the  correct  curves  and  the  altered 
curves  were  inserted  into  the  cloud-edge  computer  program  with  both 
correct  and  slightly  altered  time  information.  The  results  of  these 
sensitivity  tests  appear  in  Table  IX.  Details  of  the  sensitivity  tests 
and  complete  output  data  are  contained  in  Appendix  D. 

This  exercise  demonstrated  that  the  computer  program  is 
sensitive  to  program  Inputs.  An  error  of  more  than  about  5%  in  the 
time  between  cloud  edge  positions  or  a  systematic  error  of  about  5% 

In  the  position  of  a  cloud  edge  appears  capable  of  invalidating  program 
outputs.  In  particular,  it  is  interesting  to  note  the  similarity  between 
program  output  for  the  artificially  flattened  curve  and  that  from 
shot  #2385  (Table  V).  In  both  instances  the  computed  ejecta  velocity 
increased  with  time  and  the  angle  decreased.  (These  trends  are  the 
opposite  of  expected  behavior.  )  Also,  it  appeared  that  the  origin  was 
starting  to  move  back  in  toward  the  centerline  on  the  last  output  from 
shot  #2385,  as  it  did  consistently  with  the  altered  curves. 

Tlicse  similarities  might  Indicate  that  the  cloud  edge  data 
obtained  with  the  B&tW  300  camera  records  were  insufficiently  precis, e 
in  either  spatial  or  time  resolution  to  permit  program  solution.  On 
the  other  hauid,  the  same  result  would  occur  if  the  oarticles  were  being 
decelerated.  Particle  deceleration  computations  (Appendix  C)  indicate 
that  particles  of  10  p.  diameter  would  lose  5%  of  their  velocity  (assumed 
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TABLF  DC 

of  Cloud-Edge  Computer  Program  Sensitivity  Tests 
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Time  of  curve  A  Increased  by  5%  -Z8  to  55%  +19  to  44%  +8  to  17%  -24  to  37% 

of  (TB-TA)  moves  in  and  and  increases 

wrong  decreases 

direction 


initially  at  7  km/sec)  over  a  travel  of  about  2.4  cm.  It  doe&  not  seerri 
unreasonable  that  particles  in  the  early  ejecta  could  be  this  ;jrr,all, 
or  somewhat  smaller. 

In  order  to  reduce  the  effects  of  atmospheric  drag,  a  second 
series  of  impact  experiments  was  carried  out  at  an  ambient  oressu:* 
of  only  .09  torr.  At  this  pressure,  particles  would  have  to  approach 
0.  1  ^  in  size  before  they  would  experience  a  5%  velocity  loss  over 
6.7  cm  of  travel.  It  appears  extremely  unlikely  that  all  or  a  majority 
of  the  particles  forming  the  cloud  edge  would  be  this  small.  A  highly 
precise  image  converter  camera  system  replaced  the  framing  camera 
to  eliminate  significant  time  and  positional  errors  from  the  cloud 
photographs  in  these  experiments. 

Finally,  the  results  of  the  discrete-particle  measurements  of 
shot  #2387  can  hardly  be  ascribed  to  deceleration  caused  by  atmos¬ 
pheric  drag.  All  of  these  particles  were  in  the  range  of  0.  Z5  mm  and 
larger.  Particles  of  this  size  should  travel  more  than  20  cm  before 
experiencing  a  5%  velocity  loss  at  25.  5  torr  atmospheric  pressure. 

So,  clearly,  the  unreasonable  ejection  times  obtained  for  these 
particles  cannot  be  attributed  to  simple  particle  deceleration  due  to 
atmospheric  drag.  Although  there  was  undoubtedly  error  in  deter¬ 
mining  their  positions,  the  errors  should  have  been  random  and  not 
such  as  to  give  the  consistently  early  ejection  times  appearing  in 
Table  VI.  It  is  proposed  that  the  deceleration  of  these  large  particles 
is  caused  by  their  collision  with  smaller  particles  ejected  at  earlier 
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tliYiPB  and  themselves  decelerated.  Hence  thi  assumption  of  constant- 


velocity  travel  after  ejection  is  not  tenable. 

It  is  appropriate  here  to  review  the  assuniptlona  that  formed 
the  basis  for  the  ejecta  cloud  investigation  and  to  assess  their  validity 
in  the  light  of  information  provided  by  this  study. 

1.  It  was  assumed  that  the  angle  between  the  crater  wall  at 
the  target  surface  and  the  target  face  increases  during 
crater  growth.  The  radiographic  investigation  demon¬ 
strated  that  this  angle  remains  approximately  constant 
duriiig  almost  the  entire  crater  growth  process. 

2.  It  was  further  assumed  that  material  departs  the  crater 
along  vectors  tangent  to  the  crater  wall  at  positions 
corresponding  to  the  limits  of  crater  radius.  Measure¬ 
ments  of  cloud  base  and  throat  diameters  do  not  agree  with 
radiographic  determination  of  instantaneous  crater  radius. 
No  relationship  between  these  parameters  has  yet  been 
found. 

3.  The  ejecta  was  assumed  to  travel  in  straight  lines  at 
coixstant  velocity  after  departing  the  target  face.  It  is 
strongly  indicated  that  at  normal  range  anibient  pressures 

torr)  the  ejecta  cloud  is  significantly  decelerated  by 
atmospheric  drag  over  the  distances  required  for  the  ejecta 
cloud  8tudy%  Discrete  particle  measurements  prove  that 
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other  deceleration  forces  are  acting  on  at  least  the  late 
time  ejecta  cloud. 

4.  It  was  assumed  that  as  the  pressure  in  the  crater  dissi¬ 
pates,  the  velocity  of  ejection  decreases.  Comparison  of 
over-all  rnoven'icnt  of  the  early  ejecta  cloud  with  measured 
velocities  of  particles  in  the  late  time  cloud  grossly 
support  this  as sun'iption. 

6.  Finally,  it  was  assumed  that  at  any  instant,  particles  are 
being  ej(ictcd  at  nearly  identical  velocities.  No  supporting 
or  refuting  evidence  is  available  from  this  study. 

The  proven  Invalid:!  /  of  assumptions  1  through  3  is  felt  to 
constitute  the  reaoonc  for  the  failure  of  the  computer  program  to 
predict  crater  parameters  as  a  function  of  time.  The  assum.ptlon  of 
constant- velocity,  straight-line  ejecta  travel  is  essential  for  the 
proposed  ejecta  cloud  solution.  Therefore  the  direct  determination  of 
crater  radius  as  a  function  of  time  from  either  cloud-edge  data  or 
discrete  particle  measurements  does  not  appear  promising.  It  may  be 
possible  to  correlate  crater  radius  with  such  grose  parameters  of  the 
ejecta  clcud  as  throat  and  base  width,  but  further  investigation  will  be 
required  to  establish  a  relationship  between  these  parameters. 

Conclusions 

The  conclusions  drawn  in  the  preceding  discussion  can  be 
summarized  as  follows; 
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a.  The  crater  assumes  its  final  shape  very  early  in  the 
ejection  process,  certainly  before  3  microseconds  for  the 
experiments  being  reported.  This  shape  remalno  virtually 
unchanged  as  the  crater  continues  to  expand.  No  signi¬ 
ficant  evidence  of  "springback"  or  reduction  in  crater  size 
at  late  time  was  apparent  in  this  investigation. 

b.  The  tangent  to  the  crater  wall  at  the  target  surface  was 
found  to  make  an  angle  of  77+4  degrees  with  the  target  face 
throughout  the  course  of  crater  formation.  Hence  the 
crater  is  not  actually  hemispherical  but  is  probably  better 
represented  by  a  paraboloid. 

c.  Radiographs  appear  to  offer  the  best  means  of  observing 
crater  growth  as  a  function  of  time.  The  determination  of 
transient  crater  parameters  by  the  direct  comparison  of 
radiographs  taken  during  and  after  the  event  under  identical 
geometric  relationships  is  both  f^sasible  and  preferable  to 
any  known  absolute  technique. 

d.  Of  the  two  types  of  radiographs  employed  in  the  study, 
profile  views  of  tV>e  crater  seem  far  superior  to  plan 
views.  The  profile  view  permits  shnultaneous  determina¬ 
tion  of  both  diameter  and  depth,  while  the  plan  view  gives 
diameter  inf orrriation  only.  Moreover,  the  crater  lip 
completely  ringing  the  crater  in  the  plan  viev/  makes 
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determination  of  true  crater  diameter  more  difficult  than 


with  the  profile- view  radiographs. 

«.  For  spherical  aluminum  projectiles  impacting  thick  1100-0 
aluminum  targets  at  about  7  km/sec,  the  crater  appears  to 
follow  an  exponential  grow^’h  law  of  the  form  R  =  1-e”^*, 
where  R  is  the  ratio  of  crater  transient  dimension  to  final 
dimension.  With  t  measured  in  microseconds,  the  value 
of  k  determined  in  this  investigation  was  0.  217  1  .  003. 

f.  Use  of  three  optically  aligned  Image  converter  cameras  to 
obtain  optical  data  of  high  precision  is  entirely  feasible. 
Used  in  conjunction  with  a  Dynafax  or  similar  streaking 
camera  to  record  impact  flash  and  camera  shutter  opera¬ 
tions,  this  system  is  capable  of  providing  time  resolution 
on  tile  order  of  40  nanoseconds. 

g.  There  may  be  a  relationship  between  gross  ejecta  cloud 
parameters  (such  as  width  of  throat  and  base)  and  crater 
radius.  But  additional  evidence  vdll  be  required  to  deter¬ 
mine  what  relationship  (if  any)  exists. 

h.  The  simple  model  proposed  at  the  outset  of  the  study-- 
ejecta  particles  departing  the  crater  wall  on  a  tangent  and 

in  straight;  lines  at  constant  velocity--does  not 
agree  with  the  physical  reality.  Particle  interaction  after 
ejection  setme  quite  probable,  with  resultant  changes  in 
both  particle  speed  and  direction. 
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Recommendations 


The  following  recommendations  are  made  for  additional  work 
in  areas  investigated  in  this  study; 

a.  Additional  radiographs  should  be  made  to  observe  crater 
growth,  especially  in  the  very  early  (  <C  3  ^  sec)  and  middle 
(20-40  y.  sec)  stages.  If  the  additional  data  points  conform 
to  the  proposed  growth  law,  comparison  should  be  possible 
with  theoretical  blast-wave  and  hydrodynamic  treatments. 

b.  Further  study  should  be  made  of  the  precise  shape  of  the 
crater  during  its  formative  stage.  Again,  these  results 
should  permit  comparison  with  theoretical  predictions. 

c.  Three  image  converter  cameras  should  be  used  to  photo¬ 
graph  discrete-particle  ejecta.  In  order  to  be  able  to 
identify  and  follow  the  particles,  the  pictures  should  be 
taken  2  to  3  ^  sec  apart.  Time  and  spatial  resolution 
should  t>e  sufficient  to  determine  whether  there  is  any 
deviation  from  straight  line  travel  and  any  deceleration 
above  that  to  be  expected  from  atmospheric  drag.  These 
shots  should  he  accomplished  in  as  hard  a  vacuum  as  can 
be  maintained. 
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Appendix  A 


Computer  Program  for  Determining  Crater  Parametera 
F rom  Ejecta  Cloud  Edge  Positions 


Since  the  logic  of  the  program  is  discussed  in  Chapter  II,  this 
Appendix  will  be  limited  to  a  discussion  of  program  inputs  and  outputs. 
Because  all  measurements  were  obtained  in  English  units,  the  program 
was  written  with  inputs  and  outputs  expressed  in  these  units.  Con¬ 
version  to  the  mks-egs  system  was  made  prior  to  presentation  of  the 
data  in  the  main  body  of  this  report.  Program  inputs  and  outputs  are 
listed  below,  followed  by  a  machine  printout  of  the  entire  program. 


Program  Inputs 


X(l,  1) 


START 


XBLIM 


XCLIM 


-  Distance  {in  inches)  from  target  face  to 
starting  point  on  first  curve  (curve  A). 

-  Distance  (in  inches)  on  target  face  from  cloud 
centerline  to  assumed  point  of  origin  of 
particle  passing  through  X(l,  1). 

-  Limiting  distance  (in  inches)  from  target  face 
along  curve  E  of  acceptable  solutions.  (Out¬ 
side  the  limits  of  measurement,  the  fitted 
curves  do  not  necessarily  represent  actual 
cloud  edge  positions. ) 

-  Similar  limiting  distance  on  curve  C. 
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XBASE 

TA 

TB 

TC 

C{1.J) 

JLIM 

KLIM,  LLIM, 
ILIM 

DELTAX 

TO  LINT 

DISTOL 

ADJFAC 


-  ABsumed  tf-ne  after  impact  (in  nanoseconds) 
of  ejection  of  particle  passing  through  X(l,  1). 

-  Time  after  impact  (in  nanoseconds)  of  curve  A. 

-  Time  after  impact  (in  nanoseconds)  of  curve  B. 

-  Time  after  impact  (in  nanoseconds)  of  curve  C. 

-  Coefficients  of  curves  fit  to  cloud  edge  data. 
1(1-3)  identifies  curves  A-C.  J(l-ll) 
identifies  coefficient  of  terms  x®  through 

-  Number  of  straight-line  trajectory  outputs 
desired. 

-  Internal  iteration-number  limits. 

-  Increment  of  x-movement  desired  along  curve 
A  between  separate  straight-line  trajectory 
solutions. 

-  Acceptable  tolerance  on  iteration  routine  to 
establish  intersections  of  straight-line  tra¬ 
jectories  and  cloud-edge  curves. 

-  Acceptable  tolerance  on  iteration  routine  to 
adjust  distance -time  ratios  along  straight- 
line  trajectories  between  curves. 

-  Sensitivity  adjustment  factor  to  control 
distance  assumed  origin  is  moved  to  correct 
above  ratios. 
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Program  Outputs 


TIME 


ORIGIN 


ANGLE 

VELOCITY 

XA,  YA 

XB,  YB 

XC,  YC 


-  Time  after  impact  of  ejection  of  particle 
(nanoseconds). 

-  Lateral  (or  y-)  distance  (in  inches)  along 
target  face  from  cloud  centerline  to  point  of 
ejection  of  particle  (In  the  assumed  model, 
this  corresponded  to  crater  radius.) 

-  Angle  of  ejection  (in  degrees)  of  particle 
measured  from  plane  of  target  face. 

-  Velocity  of  ejected  particle  (in  feet/sec). 

-  Coordinates  (in  inches)  of  point  on  curve  A 
through  which  particle  passed. 

-  Same  information  for  curve  B. 

-  Same  information  for  curve  C. 
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COMPUTE  DISTANCES  BETWEEN  INTERSECTIONS 
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SHSRO'JTIME  TO  FIND  T  NTF  RE'^Ot  i  ON  OF  O'^VF  ANO  ctpatt-mt 


Appendix  B 


Raw  Data  and  Sample  Computations 


1.  Time  Computations  (Shot  2464) 

The  computations  carried  out  to  determine  times  after  im¬ 
pact  of  image  converter  camera  shutter  operations  are  presented  as 
an  example  of  similar  time- deter nrunations  made  on  all  experiments. 
Three  separate  readings  were  obtained  for  each  distance  on  both  the 
Fastax  and  Dynafax  film  records.  In  the  interest  of  brevity,  only  one 
set  of  readings  and  computations  will  be  presented.  T^  represents 
the  winker  fired  by  the  PM  tube  output, 

Fastax  Record: 


Average  film  speed  (mm/millisec)  60.  052 

Distance:  T^-Shutter  Signal  1  (mm)  .480 

T^-Shutter  Signal  2  (mm)  .  696 

TQ-Shutter  Signal  3  (mm)  .994 

Distances  were  then  divided  by  speeds  to  obtain  times: 

Time:  T^-Shutter  Signal  1  (psec)  7.993 

Tp-Shutter  Signal  2  (psec)  11.  590 

T^-Shutter  Signal  3  (psec)  16.552 

Dynafax  Record: 

Distance:  Tg-Shuttcr  Signal  1  (in)  .  0781 

Tp-Shutter  Signal  2  (in)  .1115 

TQ-Shutter  Signal  3  (in)  .1572 

Impact  -  Tq  (in)  .  0040 


Diitances  were  divided  by  (Faatax)  timea  to  obtain  film  apeed  valuea; 

Apparent  film  apeed:  TQ-Shutter  Signal  1  (in/paec)  .  0097710 

T^-Shutter  Signal  2  (in/^uaec)  .  0096204 
TQ-Shutter  Signal  3  (in^aec)  .  009497 3 
Average  film  apeed  (in^aec)  .0096296 

This  value  of  Oynafax  film  apeed  was  then  applied  to  determine  rela¬ 
tive  times; 

Impact  -  Shutter  Signal  1  8.  526  paec 

Impact  -  Shutter  Signal  2  11.994pBec 

Impact  -  Shutter  Signal  3  l6.740  )iBec 

The  time  valuea  from  each  of  the  three  separate  computations  were 

then  averaged  to  produce  the  final  valuea  of  8.  52,  12.10,  and  16.93 

paec. 

2.  Optical  Data 

a.  Discrete  Particle  Computations.  Twenty  discrete  parti¬ 
cles  were  identified  in  B  &  W  300  camera  frame  36  of  shot  #2  387  and 
followed  through  frame  44.  Their  x-  and  y-travel  over  this  time- 
span  (8.  62 ^tec)  were  determined.  Distances  divided  by  time  yielded 
X-  and  y-  components  of  velocity.  Total  distance  traveled  and  total 
velocity  were  then  determined.  The  ejection  angle  was  determined 
as  tan  ^  (Ax/Ay).  The  x-distance  of  each  particle  from  the  target 
face  divided  by  its  x-velocity  gave  the  time  of  flight.  This  value  sub- 
stracted  from  the  time  of  frame  44  (42.  2  36  ^sec  after  impact)  yielded 
the  time  of  ejection.  The  time  of  flight  of  each  particle  multiplied  by 
its  y-velocity  gave  ita  y-travel.  This  value  subtracted  from  its  y- 
distance  from  cloud  centerline  gave  its  radial  distance  of  ejection 
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(origin).  The  data  for  the  20  identified  particles  appear  in  Table  X. 
Measurements  were  made  in  English  units  and  these  are  retained  in 
the  table.  Conversion  to  the  mks  system  was  made  as  the  last  step 
prior  to  presenting  the  results  in  the  main  body  of  the  report. 

b.  Image  Converter  Camera  Data.  Coordinates  of  cloud- edge 
positions  were  read  on  the  micro'/iewer  and  converted  to  true-space 
coordinates  by  the  application  of  a  demagnification  factor  determined 
from  the  appropriate  grid  picture  (see  Fig.  34).  Film  and  real-space 
X-  and  y-coordinates  for  the  edges  of  the  outermost  or  "C"  curve  are 
presented  in  Table  XI,  The  demagnification  was  3.  34.  Distances  are 
in  inches  for  uotn  film  and  real  space. 

3.  Radiographs 

The  radiographs  from  shots  *2384,  2385  and  2387  appear  as 
Figs,  37-42. 

4,  Crater  Depth  Measurements 

Maximum  crater  depths  were  measured  from  the  plane  of  the 
undisturbed  target  face.  Two  diameter  measurements  were  made  of 
each  crater  along  perpendicular  traverses.  The  point  was  removed 
from  the  depth  gage  in  making  these  measurements  to  preclude  get¬ 
ting  an  erroneous  measurement  which  would  result  if  the  side  of  the 
point  contacted  the  crater  lip  before  the  point.  To  compensate  for  the 
absence  of  the  point,  the  shaft  diameter  was  added  to  the  measured 
value  of  crater  diameter  since,  at  the  extremity  of  each  traverse,  the 
shaft  contacted  the  crater  wall  a  distance  of  one  shaft  radius  prior  to 
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TABLE 


TABLE  XI 


Film  and  Real-Space  Coordinates  of  Cloud  Edge  "C" 
from  Inage  Converter  Photograph,  Shot  ^2464 


Film  Coordinates 

Real 

-Space  Coordinates 

Lx.  Curve  Rt 

.  Curve 

Lt.  Curve 

Rt.  Curve 

X 

Y 

Y 

X 

Y 

Y 

.0291 

-. 1332 

.  1284 

,097 

-.445 

.429 

.0541 

1353 

.1330 

.181 

-.452 

,445 

.0791 

-.1386 

.  1387 

.264 

- .  464 

.464 

.1041 

- .1426 

.1436 

,  348 

-.477 

.480 

.  1291 

-.1506 

.1504 

:  ;2 

-.504 

.503 

.:'.54i 

-.1580 

.1593 

.  r.is 

-.  526 

.533 

.1791 

-.1639 

.1667 

.599 

-.548 

.558 

,2C4l 

-.1701 

.  1740 

.683 

-.569 

.582 

.  .2291 

-.1763 

.  1814 

,  766 

-.590 

.607 

.2541 

-.1821 

.191 1 

.  850 

-.609 

.639 

.2791 

-.1903 

.20 ’.8 

.933 

-.636 

.675 

.  3041 

-.2016 

.2097 

1.017 

-.674 

.701 

.3291 

-.2147 

.2244 

1.101 

-.718 

.750 

.  3541 

-.2243 

.2335 

1,184 

-.750 

.781 

,3791 

-.2357 

.2453 

1.268 

-.788 

.020 

.4041 

-.2548 

.2549 

1.352 

-.852 

.852 

.4291 

-.2670 

.2660 

1.435 

-.893 

.890 

.4541 

-.2750 

.2780 

1.519 

-.920 

,930 

.4791 

-.2923 

.2858 

1.602 

-.978 

.956 

.5041 

-.3166 

.  3007 

1.686 

-1.059 

1.032 

,5291 

-.3361 

.3187 

1.770 

-1. 124 

1.066 

.5541 

-.3538 

.3499 

1.853 

-1.183 

1.170 

.5791 

- . 3962 

.3774 

1.937 

-1,325 

1,262 

.6041 

♦ 

.4153 

2.020 

- 

1 . 389 

*6291 

*■ 

.4467 

2.104 

- 

1.494 

*Cloud  edge  too  indistinct  to  permit  determination  of 
coordinate. 
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Fig.  37,  Radiogvaphs  of  Crater  Profile  Jluring  and 
After  Shot  #2384  (Vertical  X~ray  Setup) 
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During  Shot 


After  Shot 


During  Shot 


After  Shot 


2.  99  p  sec 


k _ J 


3.21  ;  I  s  e  r 


Fig.  40.  Radiographs  of  Crater  Diarr.eter  During  and 
After  Shot  #2385  (Horizontal  X-ray  setup) 


Fig.  41.  Radiographs  of  Crater  Profile  During  and 
After  Shot  «2387  (Vertical  X-ray  Setup) 


During  Shot 


After  Shot 


actual  shaft  centerline  position.  All  diameter  measurements  were 
made  in  the  plane  of  the  undisturbed  target  face.  The  uepth  and  two 
diameter  measurements  for  each  of  the  six  craters  appear  in  Table 
XU.  All  measurements  are  in  inches. 


TABLE  XII 

Crater  Measurement  Raw  Data 


Shot  Number 

Depth 

Diameter  1 

Diameter  2 

Z384 

.  3554 

.  6446 

.  6566 

2385 

.  3663 

,  6580 

.  6569 

2386 

.  3635 

.  6496 

.  6454 

2387 

.  3391 

.  6311 

.  6301 

2463 

.  36'.  3 

.  6490 

.  6333 

2464 

,  3413 

.  6239 

.  6354 

no 


Appendix  C 

Atn'.os p_he ric  D rag  Effects 

The  results  of  both  cloud-edge  and  discrete  particle  investiga¬ 
tions  suggested  that  the  ejecta  was  being  decelerated.  To  test  this 
hypothesis,  relationships  between  particle  size,  atmospheric  pressure, 
and  resultant  drag  were  investigated.  Work  by  Opik  (Ref  21)  and 
MacCormack  (Ref  17)  provided  the  background  for  this  investigation. 
For  the  size  particles  under  consideration,  Opik  finds  that  the  drag 
force  can  be  represented  by 

Fj^=iKSpV^  (1) 

where  is  the  drag  force,  K  is  the  drag  coefficient,  S  is  the  area 
of  minimum  convex  envelope  of  the  body,  p  is  the  atmospheric  density, 
and  V  is  particle  velocity  (Ref  21:29-38).  Opik  further  points  out  that 
the  particles  can  be  treated  as  spheres;  hence  the  equation  can  be 
rewritten 

F^  =  K  TT  r^  pv^  (2) 

where  K,  P,  and  V  have  the  same  nteanings  as  in  tht;  previous 

equation,  and  r  is  the  radius  of  the  (assu  ned)  spherical  particle. 

As  atmospheric  pressure  (and  theri'fore  density)  is  reduced, 
the  value  of  K  changes.  At  relatively  high  atmospheric  densities,  K 
has  a  value  of  0.  5.  In  this  area,  hydrodynamic  flow  prevails,  and 
the  drag  equation  assumes  the  familiar  form 
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(3) 


F 


D 


2 

r 


At  very  low  atmospheric  densities,  free  molecular  flow  exists.  Here 
K  attains  its  maximum  value  of  1.  Values  of  K  for  both  flow  regimes 
and  for  the  transition  regime  between  are  presented  in  Table  XHI. 


TABLE  Xni 

Kinetic  Thickness  and  Drag  Coefficient 


d  =  0.  75r/X» _ 

0  1 

1  0.  75 

2  0.625 

4  0.  562 

a>  0.  5 


(From  Ref  21:38) 


Here  d  is  the  average  kinetic  thickness  of  the  aircap  formed  in  front  of 
the  particle.  The  term  required  to  compute  d  represents  the  half¬ 
energy  range  of  an  air  molecule  at  the  velocity  of  the  particle,  or  the 
range  of  travel  required  for  the  molecule  to  lose  through  collisions 
one-half  its  initial  kinetic  energy.  For  a  particle  velocity  of  7  km/ 

sec  (taken  to  be  a  reasonable  estimate  of  the  velocity  of  early-time 

-8  -2 

post-jetting-phase  particles),  =6x10  gm-cm  /p,  where  p  is 
the  atmospheric  density  measured  in  c.  g.  s.  unite  (Ref  21:81). 

To  simplify  the  calculations,  it  will  be  assumed  that  nitrogen 
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molecules  comprised  100%  of  the  atmosphere  in  the  target  tank. 
From  the  ideal-gas  equation  of  atatei  atmospheric  densities  in  the 
target  tank  at  68*^F  were  computed  for  pressures  of  25.  2  and  .  09 
torr: 

-5  3 

P25  2  "  gm/cm  ; 

P  Qg  =  !•  ^8  xlO  ^  gm/cm^ 


The  appropriate  values  of  are  then 
I  6x10  ^ 

^  25.  2  ""  3.  86  xlO'^ 

,  _  6  xlO-Q 

^  .  09  ■  1.  38  xlO""^ 


1. 55  xlO 


4.  35  xlO 


cm; 


cm 


The  deceleration  of  particles  of  characteristic  dimension  (diameter) 
of  10  10  10  10  and  10  ^  cm  will  be  investigated. 

Computer  program  sensitivity  tests  (Appendix  D)  indicated 
that  a  decrease  in  particle  velocity  of  5%  would  surely  invalidate  the 
computer  output.  Therefore,  a  5%  change  in  velocity  was  selected  as 
the  limit  of  investigation.  Since  the  value  of  X.  ’  is  only  weakly 
velocity-dependent,  X.'  and  consequently  K  will  be  considered  to  re¬ 
main  constant  over  the  velocity  regime  of  interest. 

From  Newton's  Second  Law,  Eq  (2)  can  be  rewritten 

„  1  4  3  dV 

Kifr  pV  =  --  wr  P,  ^ 
a  3  b  dt 

where  and  p^  are  the  densities  of  the  target  tank  atmosphere  and 
the  (assumed  spherical)  aluminum  particle  respectively. 
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>  ...  -3K  „2 

R.«rmng  ,  — t  =  77  X  V 

b 


Since  V  =  —  •  the  equation  can  be  rewritten 


~3XPi 

V  dt  dx  dt  dx  4rPv 


Rearranging  for  integration  over  a  5%  velocity  decreaae 


V  •trCk  4  4rPh 


Solving  for  x,  the  distance  the  particle  travels  in  experiencing  a  5% 


loss  in  velocity. 


-4Pb  r 

X-—  (ln.95)(-) 


Values  of  particle  travel  for  various-size  particles  at  the  two 
pressures  of  interest  appear  in  Table  XIV. 

These  computations  are  simplified  in  that  ablation  of  the 
particle  is  ignored.  However,  the  results  for  .  09  torr  compare 
favorably  with  results  obtained  by  MacCormack  who  dealt  with  abla¬ 
tion.  He  concluded  that  in  an  atmosphere  of  0.  2  torr,  particles  with 
radii  greater  than  8  xlO  ^  cm  would  not  experience  a  10%  velocity 
loss  from  an  original  velocity  of  10  km/sec  over  a  time>span  of  5 
microseconds  (Ref  17:13).  They  are  also  in  general  agreement  with 
Carey's  conclusion  that  a  0.  25  mm  particle  will  not  decelerate 
appreciably  over  a  distance  of  30  cm  in  a  pressure  of  25.  2  torr. 

(Ref  2:111). 


Appendix  D 


Computer  Program  Sensitivity  Teats 

Detai'.&  of  the  computer  program  sensitivity  tests  (discussed 
in  Chapter  V)  appear  in  Table  XV.  All  values  are  presented  in 
English  units  as  program  output. 
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Table  XV  (Cont. ) 

Teat  Conditions  Point  Time  %  Error  Origin  %  Error  Ejection  Angle  %  Error  Velocity  %  Error 

(/laec)  (in)  (Deg)  (fps) 
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*  Program. failed  to  find  a  solution 
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